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that  the  error  probability  decreases  with  increasing  code  length.  The  multiuser 

capabilities  of  the  uplink  and  downlink  channels  were  shown  to  differ 

significantly.  It  was  shown  that  for  K  users  at  a  given  error  rate  in  the 
2 

uplink,  K  users  can  be  supported  at  the  same  error  rate  in  the  downlink. 
Fading  effects  were  analyzed  using  the  Nakagami  m-distribution  as  a  generalized 
statistical  model  of  fluctuations  in  the  received  optical  intensity.  The 


average  probability  of  error  was  derived  for  both  correlated  and  uncorrelated 
fading  in  terms  of  the  number  of  interfering  users,  the  Nakagami  fading 
parameter  and  the  average  number  of  photons  counted  in  a  bit  interval.  The 
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CHAPTER  1 

INTRODUCTION 

1 . 1  Multiple  Access  Satellite  Communications 

Multiple  access  satellite  systems  enable  a  large  number  of  ground  stations 
to  simultaneously  transmit  their  respective  information  signals  (voice,  data, 
teletype,  television,  etc.)  through  a  satellite  repeater.  An  attractive 
feature  of  the  satellite  system  is  global  accessibility.  A  single  geostationary 
satellite  is  visible  over  approximately  one-third  of  the  earth’s  surface  and  is 
theoretically  capable  of  providing  global  coverage  to  a  large  number  of  widely 
dispersed  ground  stations  [Withers,  1977].  The  information  flow  in  a  typical 
multiple  access  satellite  communication  network  would  be  similar  to  that  shown 
in  Figure  1.  To  exploit  the  unique  properties  of  global  coverage  and  multiple 
accessibility  inherent  in  satellite  repeaters,  the  various  communication  links 
using  them  must  be  separated  from  each  other.  There  are  a  variety  of  modulation 
techniques  used  to  accomplish  user  separation.  The  most  common  are  a  basic  or 
hybrid  form  of  the  following  [Wittman,  1967]. 

1)  Frequency  division  multiple  access  (FDMA):  In  FDMA  systems,  the 
bandwidth  of  the  satellite  repeater  is  divided  into  a  number  of  nonoverlapping 
frequency  bands  which  constitute  access  channels.  Each  user  is  given  exclusive 
usage  of  an  access  channel  at  a  predetermined  frequency.  Network 
synchronization  is  not  required,  but  whenever  more  than  one  signal  is  present  at 
the  repeater,  intermodulation  noise  develops.  As  the  number  of  users  increases 


Network  configuration  of  a  satellite  multiple  access 
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this  noise  builds  up  and  ultimately  limits  the  system  performance  [Schwartz,  J., 
1966].  FDMA  essentially  allows  each  ground  station  continuous  access  to  a 
predetermined  portion  of  the  available  satellite  bandwidth. 

2)  Time  division  multiple  access  (TDMA):  In  TDMA  systems,  an  interval  of 
time  (usually  called  a  frame)  is  divided  into  a  number  of  discrete  time  slots 
(which  are  considered  access  channels).  Each  user  is  given  exclusive  usage  of 
the  satellite  bandwidth,  but  transmission  is  restricted  to  one  or  more 
predetermined  time  slots  within  a  frame.  This  technique  avoids  the  problem  of 
intermodulation  noise  since  signals  from  different  users  are  never  present  in 
the  repeater  at  any  one  time.  Consequently,  TDMA  systems  require  stringent 
network  synchronization  to  avoid  the  problem  of  channel  overlap.  TDMA  allows 
each  user  exclusive  use  of  the  satellite  bandwidth  but  only  for  predetermined 
time  slots  within  a  frame  [Puente,  1971]. 

3)  Code  division  multiple  access  (CDMA):  CDMA  systems  are  designed  so  that 

multiple  users  can  communicate  simultaneously  in  the  same  frequency  band.  The 

multiple  access  capability  is  obtained  by  encoding  the  information  signal  of  a 

particular  ground  station  with  a  unique  periodic  code  sequence.  The  most  common 

methods  of  encoding  the  information  signal  are  frequency  hopping  (FH),  direct 

sequence  (DS)  modulation  or  hybrid  combinations  of  the  two  [Dixon,  1975].  In  FH 

systems,  such  as  the  TActical  Transmission  System  (TATS)  [Drouilhet,  1969],  each 

user  imposes  a  code-selected  frequency  hopping  pattern  on  his  normal 

transmission  so  that  cochannel  interference  occurs  only  when  users  are  hopped  to 

the  3ame  frequency  at  the  same  time.  To  minimize  this  interference,  FH  systems 

20 

employ  a  large  number  of  frequencies.  The  TATS  modem,  for  example,  has  2 
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discrete  frequency  choices.  FH  systems  also  require  rapid  response  frequency 
synthesizers  which  introduce  an  extra  measure  of  complexity  to  the  system. 

1 .2  Direct  Seauence/CDMA 

Direct  sequence  modulation  is  one  of  the  easiest  schemes  to  implement 
[Schwartz,  J.,  1966],  [Dixon,  1975].  In  a  simple  implementation,  the 
information  signal  is  multiplied  directly  by  a  periodic  code  sequence  and  then 
modulated  onto  a  multiple  access  carrier  using  one  of  the  more  conventional 
techniques  such  as  AM  or  FM  carrier  modulation.  A  receiver  extracts  the  desired 
information  by  correlation  techniques.  Each  user  transmits  on  the  same  carrier 
frequency,  shares  the  same  repeater  bandwidth  but  modulates  his  information 
signal  with  a  different  code  sequence.  There  is  no  requirement  for  precise 
timing  or  internal  network  control  between  the  various  users.  However, 
interference  is  generated  from  the  unavoidable  cross-correlation  between  the 
desired  user's  code  and  the  code-modulated  information  signals  from  all  of  the 
other  active  users. 

Cross-correlation  interference  can  be  minimized  by  using  specially 
constructed  code  sequences  which  are  as  nearly  random  as  possible,  i.e.,  the 
correlation  between  the  encoding  sequence  and  any  other  code  sequence  in  the 
system  is  near  zero.  An  example  of  such  sequences  is  the  family  of  Gold 
sequences  [Gold,  1967].  If  the  cross-correlation  interference  can  be  held  to 
acceptable  levels,  multiple  users  can  communicate  simultaneously  at  the  same 
frequency.  Note  that  CDMA  does  not  allow  each  user  exclusive  use  of  tne 
satellite  repeater  but  it  does  provide  each  user  continuous,  uncontrolled  access 
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and  full  use  of  the  repeater  bandwidth. 

Both  FH  and  DS  modulating  schemes  are  spread  spectrum  techniques,  i.e.,  the 
transmitted  signal  occupies  a  bandwidth  which  is  spread  substantially  wider  than 
that  of  the  information  being  transmitted.  Therefore,  space  laser 
communication,  operating  at  frequencies  where  the  bandwidth  allocations  are 
substantially  greater,  is  an  attractive  complement  for  the  wideband  signals 
generated  in  DS  and  FH  systems.  The  technical  feasibility  of  spaceborne  laser 
communications  is  well-known  for  both  heterodyne  detection  systems  [McElroy, 
1977]  and  direct  detection  systems  [Ross,  1973] .  fact,  the  system  cited  in 
the  reference  [Ross,  1973]  has  progressed  into  system  development  with  the 
scheduled  launching  of  the  first  laser  communication-equipped  satellite  in 
December  1981.  The  capability  of  the  laser  satellite  multiaccess  system  can 
also  be  extended  to  include  systems  which  provide  service  to  geographically 
separated  cells  of  several  independent  users.  The  multiple  access  receiver  for 
such  a  system  has  been  tested  by  McDonnell  Douglas  Astronautics  Company  [0l3han, 
1979]. 


1 .3  Outline  of  the  Investigation 

In  this  study,  we  analyze  the  effects  of  cross-correlation  interference  on 
the  error  performance  of  a  direct  detection  optical  CDMA  communication  system 
which  uses  direct  sequence  modulation.  Generalized  expressions  for  the  uplink 
and  downlink  error  probabilities  are  derived  in  terms  of  the  number  of 
simultaneous  users  and  the  normalized  cross-correlation  between  user  codes. 
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In  Chapter  2,  we  define  the  multiple  access  system  model,  giving  a  general 
description  of  the  optics  for  both  the  uplink  and  downlink  channels.  The 
specific  implementation  of  the  correlation  process  realized  at  the  optical 
receivers  is  also  outlined. 

In  Chapters  3  and  4,  we  provide  the  details  necessary  to  develop  the  error 
probabilities  for  both  channels.  Approximations  to  the  average  probability  of 
error  are  given  which  are  valid  for  large  signals  and  small  values  of  the 
cross-correlation  interference.  Numerical  results  are  given  for  various  signal 
parameters  and  multiple  users. 

In  Chapter  5,  the  analysis  is  extended  to  include  the  system’s  performance 
when  the  optical  channel  induces  slow  fading  on  the  received  intensity. 

Finally  in  Chapter  6,  we  conclude  with  a  summary  of  the  pertinent  results 
and  recommendations  for  continued  study. 


t 


7 


CHAPTER  2 


MULTIPLE  ACCESS  SYSTEM  MODEL 


In  this  investigation  we  have  assumed  the  configuration  of  a  multiaccess 
satellite  repeater  as  an  illustrative  model.  The  links  between  the 
transmitter-to-satellite  and  the  satellite-to-receiver  are  assumed  to  be  direct 
detection,  optical  polarization  modulation  systems  and  are  hereafter  referred  to 
as  the  uplink  and  the  downlink  channels,  respectively.  The  basic  transmitter 
and  receiver  optics  are  shown  in  Figure  2.  Spatial  acquisition,  tracking  and 
code  synchronization  are  assumed  to  be  handled  by  other  subsystems  and  are  not 
considered  in  this  analysis. 

2. 1  Uplink  Channel  Model 

In  the  uplink  channel,  we  assume  that  we  have  K  asynchronous  users  who 
transmit  a  coded  data  signal  to  a  satellite  receiver.  As  a  necessary  first  step 
in  minimizing  interference,  we  assume  that  the  users  are  sufficiently  separated 
geographically  so  that  the  field-of-view  of  a  particular  user,  when  imaged  onto 
the  detecting  surfaces  of  the  satellite  receiver,  does  not  overlap  that  of  a 
another  user.  Hence,  the  incident  optical  fields  do  not  overlap  and  we  may 
consider  the  out,-  -  of  the  satellite  receiver  as  resulting  from  the  additive 
intensities  of  each  independent  user.  In  both  the  data  signal  and  the  code 
sequence,  a  logical  zero  is  represented  by  -1,  and  a  logical  one  is  represented 
by  +1.  The  data  signal  b^t),  from  the  i-th  user,  is  a  sequence  of  constant 


Figure  2.  Transmit  and  receive  optics.  X  and  X^  are  the  emission 
rates  of  photodetector  outputs  Y(t)  and  X(t).  ¥he  code  sequence  is  a(t) 
and  the  data  sequence  is  b(t). 
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amplitude,  positive  and  negative  pulses  of  duration  T^.  The  i-th  user  is 

assigned  an  access  code  a^(t),  which  consists  of  a  periodic  sequence  of  constant 

amplitude,  positive  and  negative  pulses  of  duration  Tq.  We  refer  to  the  time 

interval  Tq  as  a  "chip"  so  as  to  distinguish  it  from  the  "bit"  which  is  usually 

reserved  for  the  time  interval  of  the  data  signal.  The  code  consists  of  L  chips 

and  we  require  that  the  length  of  the  data  bit  be  given  by  Td  s  LTc>  Generally, 

2  5 

the  code  length  is  on  the  order  of  10  to  10  chips. 


The  input  S^t)  to  the  optical  modulator  of  the  i-th  user  is  simply  the 
data  sequence  b^t)  multiplied  by  the  access  code  a^(t), 


S1(t)  •  •1(t-i1)  b1(t-i1)  (2.1) 

where  t  ^  is  a  random  delay  which  accounts  for  the  fact  that  the  users  are 
independent  and  not  time-synchronous.  The  are  uniformly  distributed  over 
(0,Td).  The  optical  modulator,  depicted  in  Figure  2a,  switches  the  laser  output 
between  right-and  left-circular  polarization  depending  on  the  polarity  of  the 
input.  The  output  intensity  of  each  user  is  constant.  Therefore,  if  the  system 
has  fC  active  users,  the  signal  at  the  receiver  consists  of  K,  equal  intensity, 
circularly  polarized  laser  beams. 
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2.2  flaxaUnK;  Channel  Model 

The  satellite  repeater  provides  simultaneous  user  access  by  broadcasting  a 
multiplexed  output  signal  which  is  the  asynchronous  addition  of  all  of  the 
active  users  on  the  uplink.  A  particular  downlink  receiver  extracts  the  desired 
signal  from  the  satellite  output  by  correlating  the  received  signal  with  a 
synchronous  copy  of  the  appropriate  multiple  access  code.  The  transmitted 
signal  on  the  downlink  differs  considerably  from  that  on  the  uplink  in  that  the 
input  S(t)  to  the  optical  modulator  of  the  satellite  takes  the  form 


K 

S(t)  s  l  a. (t-T. )  b. (t-f  )  (2.2) 

i-1  1  11  i 

where  K  is  the  number  of  simultaneous  users  and  b^t)  is  an  estimate  of  b^ ( t ) 
which  accounts  for  probable  errors  in  signal  processing  in  the  satellite.  S(t) 
is  the  addition  of  the  coded  data  signals  from  K  Independent  users.  As  such, 
the  magnitude  and  phase  (polarity)  of  S(t)  depend  on  the  number  of  simultaneous 
users.  The  satellite  encodes  this  information  by  transmitting  a  laser  beam 


which  is  both  polarization  and  intensity  modulated.  The  polarization,  right-  or 
left-circular,  depends  on  the  polarity  of  S(t)  and  the  intensity  depends  on  the 
magnitude  of  S( t) . 
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2.3  Detection  Process 

The  optical  receivers  on  both  the  uplink  and  downlink  channels  perform  the 
same  detection  process.  The  incident  intensity  is  demodulated  by  spatially 
separating  the  polarization  components  so  that  light  intensity  of  right-circular 
polarization  strikes  one  photodetector  and  light  intensity  of  left-circular 
polarization  strikes  another.  This  is  shown  schematically  by  the  separation  of 
the  upper  and  lower  photodetectors  in  Figure  2b.  The  photon  counts  from  the 
upper  and  lower  photodetectors  are  substracted  and  the  difference  is  then 
correlated  with  a  synchronous  copy  of  the  appropriate  multiple  access  code. 
This  demodulation  process  closely  parallels  the  PCM/Polarization  Modulation 
technique  described  in  the  reference  [Pratt,  1969],  except  that  we  require 
post-correlation  with  the  multiple  access  code  to  extract  the  desired  user.  The 
receiver  makes  a  decision  based  on  whether  the  correlated  difference  is  positive 
or  negative.  Other  active  users  in  the  channel  introduce  left-  and 
right-polarization  components  which  interfere  with  the  desired  signal.  This 
interference  arises  from  the  unavoidable  cross-correlation  between  the  code  of 
the  desired  signal  and  the  code-modulated  information  signals  of  the 
interferers.  The  correlation  process  generates  interference  terms  of  the  form 

fTd 

a1(t)ai(t-T1)b1(t-ii)dt  (2.3) 

0 

where  a1  is  the  access  code  of  the  desired  signal,  the  product  a^b^  is  66® 
code-modulated  signal  of  the  i-th  interferer  and  T.  is  the  bit  interval. 
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However,  if  the  cross-correlation  between  the  desired  code  and  the 
code-modulated  data  sequences  of  the  other  users  is  uniformly  small,  this 
interference  is  negligible  when  averaged  over  the  bit  interval. 

The  optical  receivers  are  subject  to  both  thermal  and  shot  noise.  For  low 
intensity,  direct-detection  receivers,  the  shot  noise  predominates  and  the 
thermal  noise  can  be  neglected.  The  shot  noise  is  generated  by  the  incident 
optical  fields,  background  radiation  and  detector  dark  currents.  The  background 
radiation,  which  is  generally  unpolarized,  is  therefore  equally  distributed 
among  the  photodetectors.  The  average  detector  dark  current  is  considered  to  be 
the  same  in  both  detectors.  Therefore,  without  loss  of  generality,  we  consider 
the  receivers  to  be  shot  noise-limited,  with  additive  background  noise  which  is 
the  cumulative  effect  of  background  radiation  and  detector  dark  current. 

For  shot  noise-limited  detection,  the  photoelectron  emissions  from  the 
detector  surfaces  are  characterized  by  a  Poisson  process  [Mandel,  19631-  If  we 
let  I ( t )  be  the  incident  optical  intensity,  then  the  probability  of  obtaining 
n  s  1,  2,...,  photons  in  an  Interval  (t,  t+t)  is  given  by 

t+T  t+T 

P(n)s[n/hf  |  I(x)dx]n  (n!)"1  expr-n/hf  j I(x)dx]  (2.4) 
t  t 


where  n  is  the  quantum  efficiency  of  the  photodetectors,  h  is  Planck's  constant 
and  f  is  the  optical  frequency. 
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If  I(t)  Is  a  known  signal  (deterministic),  (2.4)  is  the  well-known  Poisson 
counting  model.  If  I(t)  is  unknown  (random),  (2.4)  is  conditionally  Poisson  and 
the  counting  probability  is  obtained,  at  least  in  principle,  by  taking  the 
expectation  over  the  statistics  of  I(t).  The  result  is  normally  not  Poisson 
distributed. 


i 
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CHAPTER  3 


UPLINK  CHANNEL  ANALYSIS 


3.1  Conditional 


of  Error 


To  determine  the  effects  of  interfering  users  on  the  detection  error,  we 
choose  b^Ct)  as  the  desired  data  signal  and  the  K-l  remaining  data  signals  as  the 
interfering  signals.  The  receiver  is  synchronized  to  a^t),  the  desired  code. 
Since  all  time  delays  are  relative  to  the  desired  code,  we  let  t  =  0.  The 
receiver  is  shown  in  Figure  2b.  If  we  let  Y(t)  be  the  output  of  the 
left-circular  polarization  detector  and  X(t)  be  the  output  of  the  right-circular 
polarization  detector,  then  the  receiver  output  sampled  at  t=Td  is 


R  = 


T  . 
r  « 

a  (t)[Y(t)-X(t)]dt 

0 


(3-D 


where  Y(t)  and  X(t)  are  independent  Poisson  processes. 


During  the  bit  interval,  a^t)  takes  on  the  values  +1  and  -1.  If  we  define 


a+  ={  t  :  a^t)  ■  +1 ,  0  <  t  v<  Td  } 


(3.2) 


and 
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a”  ={  t  :  a1  (t)  =  -1 ,  0  4  t  Td  } 


(3.3) 


we  can  separate  (3.1)  into  integration  over  the  positive  and  negative  intervals 
of  a1 (t) 


R  s  R+  -  R 


(3-4) 


where 


R+ 


j  Y(t)dt  + 
a 


|  X(t)dt 
a' 


(3.5) 


and 


R 


(  f 

Y(t)dt  t-j  X(t)dt 
a”  a+ 


(3.6) 


The  probability  distribution  of  R,  which  is  the  difference  of  independent 
Poisson  random  variables,  is  completely  defined  if  the  mean  values  of  R+  and 
R~  are  known  [Pratt,  1969].  But  since  the  mean  and  variance  of  Poisson  random 
variables  are  equal,  we  have  shown  in  Appendix  A  that  the  probability 
distribution  of  R  can  also  be  expressed  as 
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P0(R=n)  =  exp  [-Var(R)]  { (Var(R)-E[R] )/(Var(R)+E[R]) }_n/2 
K 

•  It  |  {  Var(R)2  -  E[R]2  }  (3-7) 

I  n  i 

where  I  (•)  is  the  modified  Bessel  function  of  order  n  and  n  =  0,1,2,...,  is  the 
n 

number  of  photons  counted . 

To  analyze  the  effects  of  interfering  users  on  the  error  probability,  we 
first  derive  the  probability  of  error  conditioned  on  the  cross-correlation 
interference  given  in  (2.3).  If  =  +1  and  b^  =  -1  are  equally  likely,  the 
conditional  probability  of  error  is  given  by 


P(E  |  £  )  = 

1/2Pr(R>0  |  £,b1  =  -1)  +  1/2Pr(R<0|  C,b1  =  +1) 

+  1./4Pr(R  =  0  |  £,b1  =  -1) 

+  1/4Pr(R  =  0  (  £,b1  =  +1 )  (3.5) 


where 


£  5  {  C2,  C3,  *  *  ’  CK  } 


(3.9) 
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is  the  normalized  cross-correlation  vector  for  which  the  i-th  element  is  given 
by 


C._  =  1/Tdj  a1 (t)  ai(t-x1)  bi(t-T1)dt  (3.10) 

0 

The  last  two  terms  in  (3-8)  account  for  equality  in  the  test,  in  which  case  a 
random  choice  is  made  with  probability  1/2  of  being  correct.  The  average 
probability  of  error  is  then  determined  by  averaging  the  conditional  probability 
of  error  over  the  probability  distribution  of  the  cross-correlation  C^. 


3.2  General  Solution 

The  general  solution  to  (3-8)  is  determined  in  a  straightforward  manner  by 
obtaining  the  mean  and  variance  of  R  in  terms  of  the  cross-correlation  C^.  The 
mathematical  details  are  worked  out  in  Appendix  B;  the  results  are 


K 

E[R  I  £]  =  b  U  +  l  U.C. 

1  1  i=2  1  1 


and 


Var(R) 


2Ub  * 


K 

l  u. 

i=l  1 


(3.11) 
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where  U1 , 

Ui 

and 

are  the  average  number 

of  photons 

collected  in  the 

bit 

interval 

due 

to 

the  desired  signal, 

interference, 

and  background  noise, 

respectively. 

By 

substituting  (3-7)  into 

(3.8) ,  we  get 

an  expression  that 

can 

be  solved  in  terms  of  the  modified  Bessel  functions  and  the  Marcum's  Q  function 
[Marcum,  I960].  In  section  B.3  of  Appendix  B  we  have  shown  that  the  general 
result  for  K  users  is 


where 


P(E  J  C  )  =  1/2[ 1  -  Q(b,a)  +  Q(d,e)] 


K 

-1/4exp[-(2U.  +  l  U.  )][I.(ab)  -I.(ed)] 

D  ^  1  U  U 


(3.12) 


K 

a  =  [2U.  +  l  U. (1-C. )] 
b  i=2  1  l 


1/2 


K 

b  =  [2U.  +  2U,  +  I  U. (1+C. )] 
b  1  i=2  i  i 


1/2 


K 


c  =  [2U.  +  2U.  +  l  U.(1-C.)] 
b  1  i=2  1  1 


1/2 


K 


d  =  [2U.  +  I  U. (1+C.)] 
b  1=2  1  1 


1/2 


(3.13) 
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and  U1 ,1^  and  are,  as  we  previously  defined,  the  average  number  of  signal, 
interference  and  background  noise  photons  collected  in  the  bit  interval  T^.  The 
Marcum's  Q-function  is  defined  by  the  series  expression 


Q(a,b)  =  exp[-(a2  +  b2)/2]  l  (a/b)n  I  (ab)  (3.1U) 

n=0  n 


For  no  interference  (single  user),  the  probability  of  error  reduces  to 


P(E)s 1/2  [  1  -  Q(b,a)  +  Q(a,b)  ] 


(3-15) 


where 


a=[2Ub] 


1/2 


and 


b=[2(Ub  ♦  V])]U2 


(3.16) 


This  result,  as  we  expected,  is  identical  to  that  given  in  the  reference  [Pratt, 
1969]  for  shot  noise  limited,  PCM/Polarization  modulation  systems. 
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3.3  Special  Case :  No  Cross-Correlation  Interference 

As  we  previously  stated,  interference  is  minimized  by  using  codes  which 
have  low  cross-correlation.  The  irreducible  limit  to  cross-correlation 
interference  can  be  evaluated  by  setting  =  0  in  the  error  expression.  For 
the  case  of  K  independent  users,  the  expression  reduces  to 

P(E)=  1/2C1  -  Q(b,a)  +  Q(a,b) ]  (3-17) 


where 


a  =  [2U  ♦  l  U  ],/2 

b  i-2  1 

and 

b  =[20.  +  2U.  ♦  ?  U.]1/2  (3.18) 

b  1  i-2  1 

By  comparing  this  result  to  (3.15)  for  the  single  user  case,  we  note  that  even 
when  we  employ  codes  in  which  £.  =  0,  we  still  encounter  significant  interference 
noise.  In  fact  each  interfering  user  introduces  an  equivalent  background  count 


of  U  /2. 
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3.4  Probability  of  Error:  Large  Signal  Approximation 

Generally,  in  the  design  process,  we  are  concerned  with  obtaining  a 
reasonable  estimate  of  the  minimum  achievable  error  rate.  We  can  derive  this 
estimate  by  assuming  £.  =  0  and  large  signal  intensities  (U^  >>1).  Under  these 
conditions,  the  background  noise  can  be  neglected  and  we  can  use  the  asymptotic 
form  of  the  Q-function  [Schwartz,  M. ,  1966]  to  approximate  the  error 
probability.  For  the  case  where 


the  error  probability  is  approximately 


P(E)  :  1/2erfc{/2U  /  (/k7T  +  )} 


(3-19) 


or  for  large  K 


P(E) 


1/2erfc(U//2KU  ) 


(3.20) 


where  erfc(«)  is  the  complementary  error  function  given  by 


erfc(x) 


2 //if 


x 

[ exp(-u^)du 


0 


(3.21) 


t 
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This  approximation  to  the  error  probability  is  useful  when  making  a 
tradeoff  between  the  maximum  number  of  active  users  and  the  required  signal 
intensity  for  a  specified  error  rate.  Compare  this  result  with  phase- 
shift  keying  (PSK),  a  modulation  technique  similar  to  the  polarization 
modulation  employed  here.  The  probability  of  error  for  synchronous  detection  of 
a  PSK  binary  sequence  in  Gaussian  noise  [Schwartz,  1980]  is 


P(E)  =  1/2erfc(A//2N  ) 


(3-22) 


where  A  is  the  signal  amplitude  and  N  is  the  total  noise  power.  By  comparison, 
we  see  that,  in  the  limit  as  the  number  of  users  and  the  signal  intensities  get 

J 

very  large,  the  shot  noise  effects  tend  toward  a  Gaussian  process  and 
accordingly,  the  probability  of  error  takes  the  same  form  a3  that  for  PSK 
signalling  in  Gaussian  noise,  with  the  signal  amplitude  similar  to  U,  the 
average  number  of  signal  photons  and  N,  the  total  noise  power  similar  to  KU,  the 
sum  of  the  variances  of  the  K  independent  users.  Notice  this  equivalent 
signal-to-noise  ratio,  is  inversely  proportional  to  K,  the  number  of  users. 

In  Table  1 ,  we  give  the  generalized  error  expression  for  the  uplink 
channel,  with  some  special  cases  provided  for  comparison. 


TABU-; 


'lo^irnd  signal 


Gardner  and  Orr  [Gardner,  1979]  and  Hanlon  and  Gardner  [Hanlon,  1979] 
obtained  good  single-term  approximations  to  the  limiting  error  probability  for 
fading,  analog  CDMA  systems  in  terms  of  the  normalized  second  moment  of  the 
cross-correlation.  Since  the  shot  noise  behavior  of  optical  detection  is 
similar  to  the  fading  channel,  we  expected  that  the  error  performance  of  optical 
systems  was  also  highly  dependent  on  the  magnitude  of  the  normalized  second 
moment.  To  test  this  conjecture,  we  tabulated  the  probability  distribution  for 
the  cross-correlation  values  for  a  random  selection  of  10,  length  31  Gold  codes. 
We  computed  the  error  probabilities  using  (3*12)  for  the  two  code  pairs  having 
the  smallest  and  largest  second  moment.  This  provided  a  best  and  worse  case  for 
comparison.  The  probability  mass  functions  and  density  for  the 
cross-correlation  values  are  derived  in  the  reference  [Hanlon,  R.  C. ,  Peterson, 
G.  D. ,  and  Gardner,  C.  S.r  1980].  The  mass  functions  and  density  for  the  two 
code  pairs  having  the  smallest  and  largest  second  moment,  respectively  are  given 
in  Figure  3.  However,  in  a  practical  multiaccess  system  a  user  will  not  have  a 
priori  knowledge  of  the  other  user  codes.  Therefore,  an  exact  analysis  for  a 
particular  pair  of  codes  is  of  limited  utility  from  a  system's  point  of  view. 
We  would  primarily  be  concerned  with  the  average  performance  over  the  entire  set 
of  code  sequences.  Therefore,  we  also  show  the  average  error  performance  over 
all  possible  length  31  Gold  code  pair  combinations.  The  results  are  shown  in 
Figure  4.  As  expected,  the  best  case  is  the  code  pair  with  the  smallest  second 
moment  and  the  worse  case  is  the  code  pair  with  the  largest  second  moment.  From 
this  preliminary  analysis,  it  appears  that  a  useful  technique  for  assigning  user 
codes  is  to  consider  the  magnitude  of  the  second  moment  of  the  various 


Figure  *t .  Probability  of  error  vs.  average  signal  photons  for  the 
uplink  channel.  Best  and  worse  case  plots  are  for  length  31  Gold  code 
pairs  having  smallest  and  largest  normalized  second  moment, 
respectively.  Average  plot  represents  average  performance  over  all 
length  31  Gold  code  pairs. 
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cross-correlations  between  users.  Assign  high  duty  cycle  or  high  priority 
users,  sequences  which  generate  small  second  moments,  and  reserve  the  code 
sequences  which  generate  high  second  moments  to  low  duty  cycle  users. 

The  limiting  error  rate,  as  a  function  of  the  number  of  active  users,  was 
obtained  by  evaluating  the  error  expression  (3-17)  for  the  ideal  case  in  which 
the  background  noise  and  cross-correlation  were  set  to  zero.  Results  for 
systems  with  up  to  10  users  are  given  in  Figure  5  for  the  uplink  channel.  This 
plot  shows  the  degrading  effects  of  increasing  the  number  of  users  on  the  uplink 
channel.  Each  additional  interferer,  as  we  previously  noted,  has  the  effect  of 
adding  an  equivalent  photon  count  equal  to  U/2  to  the  interference  noise,  even 
though  the  cross-correlation  between  users  is  negligible.  Therefore,  as  the 
number  of  users  increases,  much  larger  signal  intensities  are  required  to 
maintain  satisfactory  error  rates.  The  curve  also  gives  a  reasonable  estimate 
of  the  minimum  achievable  error  rates  for  a  given  number  of  simultaneous  users. 
It  also  gives  a  relative  indication  of  how  much  the  signal  count  has  to  be 
boosted  to  maintain  the  error  rate  goals  as  the  system  capacity  increases.  A 
significant  point  is  that  the  error  rate  does  not  saturate  at  some  fixed  value. 
It  appears  that  if  the  users  are  capable  of  continually  boosting  their 
transmitted  power,  the  error  rate  can  be  continually  lowered.  This  is  an 
advantage  the  uplink  has  over  the  downlink  because  the  satellite-to-user 
downlink  is  power-limited  to  the  maximum  output  of  the  satellite  transmitter  and 
can  not  continually  boost  its  output  power  to  achieve  lower  error  rates. 
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As  the  code  sequences  get  longer,  their  properties  more  closely  approach 
the  properties  of  the  random  codes.  Thus,  we  expect  the  longer  code  sequences, 
with  correspondingly  low  cross-correlations,  to  improve  system  performance. 
Figure  6  shows  the  effects  of  increasing  code  length  on  the  probability  of  error 
for  a  system  of  two  users.  In  this  curve  the  background  noise  is  set  to  zero  to 
highlight  the  effects  of  increasing  the  code  length.  The  plots  are  for  code 
lengths  31.  127,  and  infinite.  The  curve  for  infinite  code  lengths  corresponds 
to  setting  the  cross-correlation  to  zero,  which  is  the  case  for  purely  random 
codes.  As  the  code  length  increases,  the  error  performance  correspondingly 
improves  and  error  rates  close  to  the  irreducible  error  rate  (that  for  infinite 
code  lengths)  can  be  achieved  using  moderately  long  code  sequences.  For 
instance,  the  performance  of  length  127  code  sequences  is  very  close  to  that  for 
infinite  codes. 


UPLINK  CHANNEL  (2  USERS) 
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CHAPTER  4 

DOWNLINK  CHANNEL  ANALYSIS 

4. 1  Conditional  Probability  of  Error 

Since  the  receivers  on  the  downlink  channel  perform  basically  the  same 
demodulation  process  as  is  performed  in  the  uplink  channel,  the  conditional 
error  probability  is  obtained  from  the  error  expression  (3-8)  and  determined  by 
substituting  the  mean  and  variance  of  the  receiver  output  resulting  from  the 
downlink  modulation  input  S(t)  defined  in  (2.2).  From  the  derivation  in 
Appendix  B,  the  conditional  mean  and  the  variance  are 


K  . 

EtR  1  £]  a  b,U  +  U  F  C.  (4.1) 

1  i=2  1 


where 


(4.2) 


and 


Var(R)  =  2U.  +  K'U 
b 


(4.3) 
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where  C.  is  defined  by 
1 


C.  =  1/T.  a  (t)  a.(t-T.  )b.  (t-xjdt 
1  d  1  o  1  i  i  i  i 


(4.4) 


and 


K'  =  1 /T , 


d  K 


l  a  (t-x  )  b  (t-x 
i=l  1  1  1  1 


.)  dt 


(4.5) 


For  K  users  the  conditional  error  probability  from  Appendix  B  is 
P(E  I  C  )  =  1/2[  1  -  Q(b’,a')  +  Q(d',c')  ] 

+1/4  exp[-(2U,  +  K ' U ) ]  In(a'b’) 
b  +  0 

-1/4  exp[-(2Ufe  +  K^U)]  I0(c’d')  (4.6) 


where 


*  1  /p 

a'  =  [2U.  +  U(K '  -  K)  +  U  l  ( 1-C. ) ] 
b  +  i=2  1 

K  -  1/2 

b'  =  [2'J,  +  2U  +  U(K '  -  K)  +  II  (1+C.)] 

D  +  •  “  o  1 
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K 

=  [2U,  +  2U  +  U(K'  +  K)  +  U  7  (1-C,)] 
b  -  i=2  1 


1/2 


K 

d'  =  [2Ufe  +  U(K '  +  K)  +  U  [  (1+Ci)] 

i=2 


1/2 


(4.7) 


The  designations  K ^  and  K*  denote  the  conditional  values  of  K'  when  b ^ ( t )  =  +1 
and  b ^ ( t )  =  -1,  respectively. 


In  the  case  of  no  interfering  signals,  the  term  K'  is 


K' 


(4.8) 


K'  reduces  to  one  since  the  product  | a ^ ( t— t 1 ) b 1 ( t— ) | is  identically  one  for  all 
values  of  a1  (t-x^  5^  (t-T  ) .  The  error  probability  for  this  case  is 


P(E)  =  1/2[  1  -  Q(b,a)  +  Q(a,b)] 


(4.9) 


where  the  arguments  a  and  b  of  the  Q-function  are  the  same  as  those  given 
(3.16).  This  result,  as  we  expected,  is  identical  to  the  uplink,  single  user 


case . 
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4.2  Special  Case :  Single  Ir.terferer 

For  an  arbitrary  number  of  users,  K'  defines  a  complex  statistical 
relationship  between  the  code  sequences  a^(t)  and  the  cross-correlation  C^; 
however,  we  can  obtain  simplified  expressions  for  K’  under  some  limiting 
conditions.  For  example,  for  two  users  we  have  shown  in  Appendix  C  that  K' 
reduces  to  the  particularly  simple  form: 


K'  =  1  +  b1C2  (4.10) 

The  resulting  error  probability  is 

P(E\C2)  =  1/2[1  -  Q(b',a')  +  Q(a',c')] 

+1/4  exp{  -[2Ub  +  (1  +  C2  )U]}  I0(a'b’) 

-1/4  exp{  -t2Ufe  +  (1  -  C2  )U]}  I0(a'c')  (4.11) 


=  [  2UK  ] 


1/2 


where 
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*  t  /? 

b*  =  [2Ub  +  2U(1  +  C2  )]  c 

c*  =  [20.  +  20(1  -  C,  )]1/2  (4.12) 

D  2 


Now  to  obtain  the  irreducible  limit  to  cross-correlation  interference  effects 
for  two  users,  we  evaluate  the  expression  above  for  =0.  The  error 
probability  feduces  to  the  same  result  obtained  for  the  single  user  case: 


P(E)=1/2[1  -  Q(b,a)  -  Q(a,b)]  (4.13) 


where  the  arguments,  a  and 

t>, 

of  the 

Q-function 

are 

given  by  (3-16). 

This 

implies  that  two  users 

employing 

code  sequences 

in 

which  C2  =  0, 

can 

theoretically  communicate 

in 

the  downlink  channel 

without 

interfering  with 

each 

other.  However,  this  is  a  fortuitous  consequence  of  the  modulation  technique 
and  does  not  hold  for  more  than  two  users. 


For  =  0,  K'  is  no  longer  coupled  to  the  cross-correlation  values,  as  is 
shown  explicitly  for  the  single  user  case,  and  may  be  treated  as  a  random 
variable  which  is  simply  the  integral  over  the  absolute  value  of  the  sum  of  the 
product  terms  a  b^.  The  probability  of  error  for  this  case  is 


P(E)  =  1/2[1  -  Q(b' ,a ' )  +  Q( a'  ,b')3 


(4.14) 
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where 


a'  =  [2Ub  +  U(K '  -  1)] 


1 


b'  =  [2Ub  +  U(K'  +  1)] 


(4.15) 


4.3  Probability  of  Error :  Large  Signal 


We  can  obtain  an  estimate  of  the  minimum  achievable  error  rate  for  the 
downlink  channel  by  setting,  as  before,  C  -  0  and  assuming  large  signal 
intensities.  Under  these  assumptions  the  background  noise  is  negligible,  and  K' 
is  independent  of  the  cross-correlation.  Using  the  asymptotic  form  of  tne 
Q-function  we  can  approximate  the  error  probability  by 


PCE)  I  1/2erfc{/2u" /  C/K'  +  l  +  /FIT )) 


(4.16) 


or  for  large  K' 


P(E)  :  1/2erfc{U//2k'U  ) 


(4.17) 


37 


K'  simplifies  considerably  when  evaluated  under  the  assumption  that  the 
cross-correlation  is  zero,  because  as  we  previously  stated  for  this  limiting 
case,  we  can  treat  K'  as  an  independent  random  variable.  The  product  terms  in 
(4.5)  are  mutually  independent,  therefore,  the  general  expression  for  K'  has  the 
form  of  a  time  average  of  the  absolute  value  of  the  sum  of  finitely  many 
independent  random  variables  which  take  on  the  values  +1  and  -1.  Assuming  that 
the  system  is  stationary,  we  can  model  the  summation  in  the  integrand  of  K *  as 
resulting  from  the  chance  outcome  of  a  random  walk  process  [Papoulis,  1965],  and 
take  K'  as  the  statistical  average  of  the  absolute  value  of  the  path  length  of 
the  random  walk  after  K  random  steps  have  been  taken,  where  K  represents  the 
number  of  simultaneous  users  in  the  system.  The  details  of  this  derivation  are 
contained  in  Appendix  C;  the  results  are 


K'  =  K 


(1/2)' 


for  K  even 


|K/2 


(4.13) 


and 
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K'  =  K 


(1/2)*-1 


for  K  odd 


UK-D/2  J 


(4.19) 


For  a  large  number  of  users,  we  can  apply  Stirling's  approximation 
[Abramowitz,  1965],  to  the  binomial  coefficients  above  and  simplify  K'  to 


K’  =/2K/ * 


(4.20) 


Comparing  the  large  signal  approximation  of  the  error  probability  for  the 

uplink  channel  with  that  for  the  downlink,  we  note  that  both  systems  degrade  in 

proportion  to  the  increase  in  users;  but  since  K'  increases  in  proportion  to 
1/2  . 

K  ,  the  downlink  channel  will  invariably  perform  significantly  better; 
particularly  for  large  K. 

In  Table  2,  we  give  the  generalized  error  expression  for  the  downlink 
channel,  with  some  special  cases  provided  for  comparison. 
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4.4  Numerical  Results 

To  further  test  the  conjecture  that  the  magnitude  of  the  second  moment  of 
the  cross-correlation  values  is  a  significant  factor  in  determining  the  error 
rate  performance,  we  evaluated  (4.11)  using  the  same  length  31*  Gold  code  pairs 
that  were  used  in  the  uplink  channel  analysis.  The  probability  of  error  was 
evaluated  for  the  two  code  pairs  having  the  respective  largest  and  smallest 
second  moment.  The  results  are  shown  in  Figure  7.  The  average  error 
performance,  calculated  using  the  probability  distribution  for  the  average  over 
all  possible,  length  31  Gold  code  pair  combinations  is  also  shown  for 
comparison.  As  with  the  uplink,  the  best  case  is  the  code  pair  with  the 
smallest  second  moment  and  the  worst  case  is  the  code  pair  with  the  largest 
second  moment.  Hence,  for  this  random  selection  of  length  31  Gold  codes,  the 
dependence  of  the  error  rate  performance  on  the  magnitude  of  the  second  moment 
of  the  cross-correlation  is  consistent  for  both  the  uplink  and  downlink 
channels. 

Comparing  Figure  7,  which  gives  the  performance  of  the  downlink,  with 
Figure  4,  the  corresponding  performance  for  the  uplink,  we  note  that  the 
downlink  channel's  performance  is  significantly  better  than  that  of  the  uplink. 
Thi3  is  not  unexpected,  for  as  we  showed  earlier  in  the  asymptotic  error 

expressions,  the  equivalent  SNR  in  the  uplink  is  asymptotically  proportional  to 

1  /2 

1/K,  while  in  the  downlink  it  is  asymptotically  proportional  to  1/  K  .  The 
performance  comparison  in  Figure  8  illustrates  the  wide  disparity  in  the  uplink 
and  downlink  channel  performance.  From  a  system's  point  of  view,  this  suggests 
that  for  every  K  users  on  the  uplink,  the  downlink  has  the  capability  of 
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Figure  7.  Probability  of  error  vs.  average  signal  photons  for  the 
downlink  channel.  Best  and  worse  case  plots  are  for  length  31  Gold  code 
pairs  having  smallest  and  largest  normalized  second  moment, 
respectivel y .  Average  plot  represents  average  performance  over  all 
length  'jl  Gold  code  pairs. 
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Data  are  given  for  systems  of  two  users  employing  uncorrelated  code 


M3 

2 

supporting  K  users  at  a  comparable  error  rate.  This  of  course  assumes  that  we 
are  operating  under  low  error  rate  conditions  in  which  we  can  neglect  the 
processing  errors  onboard  the  satellite. 

The  limiting  error  rate,  as  a  function  of  the  number  of  active  users,  was 

obtained  for  the  downlink  by  evaluating  (4.14)  for  the  ideal  case  in  which  the 

background  noise  and  cross-correlation  were  set  to  zero.  Results  for  systems 

with  up  to  20  users  are  given  in  Figure  9  for  the  downlink  channel.  This  plot 

shows  the  degrading  effects  of  increasing  the  number  of  users  in  the  downlink 

channel.  However,  in  comparison  to  the  uplink,  the  downlink  performs 

substantially  better.  The  curve  shows  that,  for  comparable  photon  counts  and 

error  rates,  the  downlink  channel  can  support  substantially  more  users  than  the 

uplink.  At  very  large  intensities,  we  have  shown  previously  that  for  every  X 

o 

users  on  the  uplink,  the  downlink  has  the  capability  of  supporting  X  users  at  a 
comparable  error  rate. 

We  also  evaluated  the  downlink  error  performance  for  length  31  and  127  code 
sequences,  to  compare  with  the  corresponding  results  for  the  uplink  channel.  We 
used  the  3ame  probability  densities  for  the  cross-correlation  values  that  were 
used  in  the  uplink  analysis.  The  results  are  shown  in  Figure  10.  The  plot  for 
an  infinite  code  length  corresponds  to  setting  the  cross-correlation  to  zero, 
which  is  a  property  of  the  purely  random  codes.  As  shown  in  the  figure,  the 
error  performance  improves  as  the  code  length  increases.  This  behavior  is 
consistent  with  that  for  the  uplink,  except  at  smaller  error  rates.  Mote  also 
that  error  performance  close  to  that  for  the  random  codes  can  be  achieved  by 
using  moderately  long  codes. 


AVERAGE  SIGNAL  PHOTONS  (U) 
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CHAPTER  5 

SYSTEM  ANALYSIS  FOR  THE  FADING  CHANNEL 

5. 1  Model  for  the  Fading  Channel 

An  optical  signal  propagating  through  the  earth's  atmosphere  can  be 
severely  degraded  by  atmospheric  turbulence  [Lawrence,  1970].  Thi3  turbulence 
is  caused  by  random  variations  in  the  refractive  index  of  the  atmosphere.  The 
interaction  of  the  optical  signal  with  the  turbulent  medium  leads  to  random 
amplitude  and  phase  variations  in  the  optical  signal  [Kerr,  1970].  These 
effects  produce  a  multipath  phenomenon  that  gives  rise  to  random  fading  of  the 
received  optical  intensity.  Theoretical  and  experimental  studies  of  wave 
propagation  in  the  atmosphere  have  shown  that  fading  at  optical  frequencies  is 
governed  by  the  lognormal  probability  density  [Tatarski,  1961].  The  strength  of 
fading  is  defined  in  terms  of  the  log- intensity  variance,  which  increases  as  the 
strength  of  fading  and  length  of  the  optical  path  increase.  Measured  data 
[Tatarski,  1961]  on  terrestrial  optical  paths  yield  values  of  the  log-intensity 
variance  ranging  from  .042  for  a  2000  meter  path  to  .0073  for  a  250  meter  path. 
These  results  are  typical  of  earth-to-satellite  links  where  the  primary 
disturbing  medium  is  the  first  few  kilometers  of  the  lower  atmosphere  [Kennedy, 
1969].  Optical  fading  generally  occurs  at  rates  less  than  500Hz  [Churnside, 
1980] . 


47 


If  the  optical  field  incident  on  the  photodetectors  in  an  optical  receiver 
fades  randomly  in  time,  the  photodetector  output  is  described  as  a  conditional 
or  doubly  stochastic  Poisson  process  [Snyder,  19751.  In  other  words,  the 
photodetector  output  is  Poisson  distributed  but  with  an  intensity  which  is 
itself  a  random  process.  The  counting  statistics  for  a  conditional  Poisson 
process  can  be  obtained  by  taking  the  expectation  of  the  conditional  Poisson 
count  with  respect  to  the  statistics  of  the  optical  intensity.  For  example,  in 
a  stationary  system,  the  probability  of  obtaining  n  =  0,  1,2,  .  .  .  ,  photon 
counts  in  an  interval  (t,  t  +  t)  is 


t+x 


t+T 


P(n)  =  Ex[(n!)_1  { 


A  (x)dx}  exp  {-J  A(x)dx)}] 
t  t 


(5.1) 


where  E  [ *  1  is  the  expectation  with  respect  to  A. A  (t)  is  given  by  the  formula 
[Pratt ,  19691 , 


A(t)  =  nKt)/hf 


(5.2) 


where  I ( t )  is  the  randomly  fading  optical  intensity,  h  is  Planck's  constant,  f 
is  the  optical  frequency  and  n  ,  (n  <  1  )  is  the  photodetector  quantum 
efficiency . 
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A  well-known  statistical  model  of  fading  signals  is  the  Nakagami 

m-distribution  [Nakagami,  i960].  If  the  magnitude  of  a  fading  optical  field  is 

m-distributed ,  the  distribution  of  the  intensity  variations  is  the  familiar 
2 

gamma  or  chi,  distribution.  This  distribution  is  used  as  a  statistical  model 

£tQ 

of  speckle  noise  in  laser  systems  [Gardner,  1977],  [Goodman,  1975].  The  photon 
counting  probability  corresponding  to  the  gamma  distributed  intensity  is  the 
negative  binominal  distribution.  The  negative  binominal  distribution  includes 
as  special  cases,  the  Bose-Einstein  distribution,  m  =  1  and  the  Poisson 

distribution,  m  -+■  ».  Under  conditions  when  the  gamma  distribution  is  used  to 
approximate  the  lognormal  distribution  [Yoshikawa,  1973],  the  negative  binominal 
is  also  a  good  approximation  to  the  counting  statistics  for  a  lognormally 
distributed  intensity.  In  this  analysis  we  use  the  Nakagami  m-distribution  as  a 
generalized  statistical  model  of  the  fading  channel.  We  assume  that  fading  is 
due  solely  to  the  properties  of  the  channel,  as  opposed  to  the  nature  of  the 
optical  source,  i.e.,  the  optical  field  originates  from  a  highly  coherent  source 
but  is  perturbed  by  random  variations  in  the  channel  properties.  We  also  assume 
that  the  background  noise  and  cross-correlation  interference  are  negligible  so 
as  to  facilitate  the  mathematical  analysis. 


5.2  Probability  Density  for  the  Optical  Intensity 

If  the  magnitude  of  the  optical  field  is  Nakr.gami  m-distributed,  the 
optical  intensity  (magnitude-squared)  is  gamma  distributed  with  probability 
density  function  [Goodman,  1975] 


1 


49 


p( I)  =  (m  /<I>)m{Im-1  /  F(m) }  exp  {-m(I/<I>)} 


(5.3) 


where  I  is  the  instantaneous  optical  intensity,  <I>  is  the  average  intensity  and 
the  factor  m  is  defined  as  the  inverse  normalized  variance  of  I, 


m  =  <I>2/<  (I  -  <I>)2  > 


(5.4) 


The  m-distribution  is  completely  specified  by  the  two  parameters  <I>  and  m. 
For  a  given  average  intensity  <I>,  m  can  be  selected  to  describe  the  strength  of 
fading.  The  strength  of  fading  is  inversely  proportional  to  m.  The  photon 
counting  probability  corresponding  to  the  gamma  distributed  optical  intensity  is 
the  negative  binominal  distribution  [Saleh,  1979], 


P(10  =  F  (m+k)/k'.r(m)  [1  +  <I>/m]~m[  l  +  m/<I>  ]~k 


(5.5) 


The  mean  and  variance  of  (5-5)  are  given  by 


E[k]  =  <I> 

and  (5.6) 

Var(k)  =  <I>  +  <I>2/m 


-1 

J 
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In  the  limit  as  m-*00  , 

T (m+k)/r(m)  :  mk 

and  (5.7) 

[1  +  <I>/m]-m  ~  exp(-<I>) 

Substituting  these  approximations  into  the  limiting  form  (m-«°  )  of  (5.5),  we 
get 


POO  =  exp( -<I>)  <I>k/k! 


(5.8) 


which  is  the  well-known  Poisson  counting  model  for  the  no  fading  case  analyzed 
in  Chapter  3.  For  the  special  case,  m  =  1,  the  negative  binominal  distribution 
reduces  to  the  Bose-Einstein  distribution 


P(k)  =  [1  +  <I>]-1[<I>/(  1  +  <I>  )]k 


(5.9) 


which  is  a  statistical  model  of  a  narrow  band  Gaussian  optical  field  [Saleh, 


1978]. 
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5.3  Unlink  Channel  Analysis 

5.3.1  Correlated  fading 

It  appears  that  the  simplest  approach  to  deriving  the  average  probability 
of  error  for  the  fading  case  is  to  generalize  the  error  expressions  derived  in 
Chapter  3  to  probability  of  error  expressions  conditioned  on  the  fading 
intensities  and  then  remove  the  conditioning  by  averaging  over  the  fading 
statistics  of  the  individual  users.  The  average  probability  of  error  is 
obtained  by  evaluating 


P(E)  = 


• 

c 

, 

P(u1fu2 


V 


P(E|U1,U2,  .  .  .  Ur)  dU1dU2  .  .  .  dUK 


(5.10) 


where  P ( E ]  U 1 , U2 , .  .  .  UR)  is  the  conditional  probability  of  error  and  pCU^Ug,. 

.  UR)  is  the  joint  density  for  U^U^.  .  .  U^,  the  conditional  photon  counts 
of  the  K  system  users. 

However,  if  we  have  equal  power  users,  the  background  noise  and 


cross-correlation  interference  are  neglected  and  the  fading  is  correlated 
between  users,  the  average  probability  of  error  is  obtained  with  a  single 
integration  over  the  probability  density  for  the  optical  intensity 
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P(E)  =  p(U)P(£IU)du 


(5.11) 


where, in  this  case,  we  set 

Va2 . UKs0 

and  use  the  error  expression  derived  in  Chapter  3  for  P(EHJ).  The  probability 
density  for  the  optical  intensity  is  p(U). 

However  for  uncorrelated  fading,  a  condition  under  which  each  user  fades 
independently,  it  is  difficult  to  derive  a  useful  form  of  the  probability  of 
error  using  this  approach.  The  difficulty  arises  from  the  fact  that  this 
approach  requires  K-fold  integrations  involving  the  joint  statistics  describing 
the  random  fading.  To  circumvent  this  problem,  we  resorted  to  an  alternative 
formulation  of  the  error  expression.  First,  we  analyze  the  simple  case  of 
correlated  fading  on  the  uplink,  under  the  conditions  of  equal  power  users, 
neglecting  the  background  noise  and  cross-correlation  interference.  For  this 
case,  we  generalize  the  probability  of  error  in  (3.18)  to  a  conditional 
probability  of  error 


P(ElU)  =  1/2[  1  -  Q(b,a)  +  Q(a,b)  ] 


(5.12) 


where  the  arguments  of  the  Q-function  are 
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a  =  [  (K  -  1)U  ]1/2 


and 


b  =  [  (K  +  DU  ]1/2 


(5.13) 


The  photon  count  U  reflects  the  equal  power,  correlated  user  case,  for  which 


U 


1 


U 


The  average  probability  of  error  is  obtained  from 


P(E) 


P(E( U)  p(U)  dU 


(5.14) 


where  p(U)  is  the  gamma  probability  density  in  the  parameter  U.  In  Appendix  D 
we  give  the  mathematical  operations  required  to  simplify  the  error  expression, 
the  final  result  is 


P(E)  =  [  1  +  K<U>/m  ]"m 


•l  C(K  -  1 )/2K]P  op  [(m)  /p! ] 
p=0  p 


•  t  [(K  +  1)/2K]J  <5J 

j-o 


[(..p)/jnv3iP 


(5.15) 


where 


5H 


5  =  [  1  +  m/K<U>  ] 

m  is  the  fading  parameter,  <U>  is  the  average  number  of  photons,  K  the  number  of 

simultaneous  users  and  y .  =  1/2  for  j  =  p  and  is  1,  otherwise. 

J » P 

5.3*2  Uncorrelated  fading 


For  the  case  of  uncorrelated  fading,  we  are  confronted  with  the 
difficulties  involved  in  deconditioning  the  error  expression  with  respect  to  the 
fading  statistics  of  each  individual  user.  To  recast  the  problem  in  a  more 
convenient  form,  recall  from  the  analysis  in  Chapter  3  that  the  receiver  output 
R  (Figure  2b),  can  be  written  as  the  difference  between  two  photon  counts,  R+ 
and  R  .  If  a  +1  is  transmitted,  R+  is  the  sum  of  the  desired  signal  count  S 
(which  is  synchronized  to  the  data  bit)  and  an  interference  count  I+. 


R  =  S 


(5.16) 


I*  represents  spurious  photons  counted  on  the  upper  photodetector  due  to  the 
presence  of  interfering  users  which  are  not  synchronized  to  the  data  bit.  On 
the  other  hand,  R  is  simply  the  corresponding  count  which  appears  on  the  lower 
photodetector,  also  due  to  the  interferers 
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R  =  I 


(5.17) 


For  no  fading,  S  and  I+  are  independent  Poisson  random  variables;  therefore,  the 
sum  H+  =  5  +  I+  is  a  Poisson  random  variable.  I”  is  also  a  Poisson  random 
variable,  independent  of  S  and  I+.  Consequently,  the  receiver  output  R  is 
analyzed  as  the  difference  between  two  independent  Poisson  random  variables. 
However,  if  we  allow  each  user  to  fade  independently  and  if  the  fading  is  slow 
enough  so  that  we  may  consider  it  essentially  constant  over  the  bit  period,  then 
we  must  account  for  the  fact  that  the  interference  on  the  upper  and  lower 
photodetectors  is  correlated  by  virture  of  the  fading  channel.  As  a 
consequence,  we  cannot  treat  the  receiver  output  as  the  difference  of 
independent  Poisson  random  variables. 

To  derive  an  alternative  expression  we  note  that  since  it  is  equally  likely 
that  b  =  +.  1 ,  the  average  probability  of  error,  defined  at  the  receiver  output 
R,  can  be  written 


P(E)  =  Pr( R+<  R")  +  1/2  Pr(R+  =  R-) 


(5.13) 


which,  using  (5.16)  and  (5.17),  can  also  be  expressed  in  terms  of  the  desired 
signal  and  interference 


P(S)  =  Pr(S  <  Z)  +  1/2Pr(S  =  Z) 


(5.19) 
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where  Z,  given  by 

Z  =  1“  -  I+  (5.20) 

is  the  difference  between  the  interference  counts  on  the  upper  and  lower 
photodetectors.  The  last  terms  in  the  error  expressions  above  account  for 
equality  in  the  decision  test,  whereby  we  make  a  random  decision  with  a 
probability  of  1/2  being  incorrect. 

In  terms  of  counting  probabilities,  the  average  probability  of  error  is 

<»  k 

P(E)  =  l  P2(k)  l  Ps(J)  Yj(k  (5.21) 

k=-oo  j=0 

where  Pg  and  Pz  are  the  counting  probabilities  for  the  desired  signal  count  S 
and  the  interference  count  Z. 

The  factor 


otherwise 


(5.22) 
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accounts  for  equality  in  the  decision  test.  At  this  point  we  must  determine  the 
unconditioned  counting  statistics  for  the  desired  signal  and  the  interference. 

To  determine  Ps  we  note  that,  since  the  system  is  synchronized  to  the  code 
of  the  desired  signal,  S  is  Poisson  distributed  given  the  fading  optical  field 

Td  Td 

Pa|^(j)  =  (J!)'1  (I  A  (x)dx) ^  exp  [-f  A(x)dx]  (5.23) 

0  '  0 

where  the  counting  rate  \(x)  is  the  random  process  proportional  to  the  intensity 

and  T,  is  the  length  of  the  data  bit  interval.  If  the  fading  rate  is  slow  when 
d 

compared  to  the  data  bit  interval,  as  will  be  the  case  for  data  rates  >  500  Hz, 
we  can  assume  A(x)  is  approximately  constant  over  the  integration  period  and 
remove  it  from  within  the  integral.  If  we  let 


f  d 

Ut  =  j  A(x)dx  I  \Td 

0 


(5.24) 


then,  U1  is  a  random  variable  which  reflects  the  fluctuating  behavior  in  the 
optical  intensity  of  the  desired  signal.  The  counting  probability  Pg  is  then 
determined  by  integrating  the  conditional  Poisson  count  over  the  gamma 
probability  density  pUJ^).  From  (5.5)  this  gives 

Pg(J)  =  "( j+tn)/  j !  T(m)  [  1  +  <U1>/m  ]"“  [  1  +  m/<U1>]“j 


for  j  >.  0  . 


(5.25) 
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where  <u1>  is  the  average  number  of  photons/bit  and  m  is  the  Nakagami  fading 
narameter . 

The  interferers,  which  are  not  synchronized  to  the  desired  signal,  generate 
spurious  photon  counts  on  both  the  upper  and  lower  photcdetectors.  The 
interference  count,  conditioned  on  both  the  fading  statistics  and  the 
cross-correlation  values  for  each  user,  is  distributed  as  the  difference  of  the 
Poisscn  random  variables,  I+  and  I-.  The  conditional  probability  distribution  of 
the  interference  count  is  [Pratt,  1969] 

P,,„  (k)  =  exp  {- (U  +  U  )}  (U  /  U  )'k/2 

z\u,c  +  -  + 

I|k  j(2vu+u_)  (5.26) 


where 

II  =  {  U2,  u3,  .  .  .  uK) 

are  the  random  variables  defining  the  average  photon  count  for  the  interferers 
and 

C  =  {  C2,  Cy  .  .  .  CK) 

are  the  cross-correlations  between  the  desired  signal  and  the  interferers.  U+ 
and  U  are  the  conditional  average  values  of  the  total  interference  counts  on 
I+  and  l”,  respectively. 
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U+  =  EC  I+  |  C.i!] 

and 

U_  =  E[  I"  |  £,ii]  (5.27) 

To  define  P  .  ,  we  must  determine  the  conditional  values  of  U  and  U  .  This 

z|u_,c_  + 

can  be  determined  from  a  simple  analysis  of  the  receiver  output  for  the  no 
fading  case  i.e.,  the  case  for  which  ji  is  deterministic.  For  K  users,  in  the 
absence  of  background  noise,  we  can  express  (3. 11)  as  the  conditional  mean  and 
variance 


K 


:[h|u,c]  =  b  u  ♦  V  c.  u. 

1  1  i-o  1  1 


and 


K 

Var(Rlil)  =  '  U.  (5.29) 

i=l  1 

where  b^  is  the  data  bit  of  the  desired  signal  and  U ^  is  the  corresponding 
average  photon  count.  is  the  average  photon  count  of  the  i-th  int erferer  and 
is  the  cross-correlation  between  the  desired  signal  and  the  i-th  interferer. 

From  results  in  Appendix  A  we  can  write 
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E[R+|C,H  ]  =  1  /2 {  Var( R  [u)  +  E[r|c,UJ  } 


and 


(5.29) 


E[R“|i,I  ]  =  1/2{  Var(R  lU)  -  E[R|C,n]  } 


By  substituting  (5.28),  evaluated  for  b  =  +1 ,  into  (5.29),  we  get 


K 

E[R+  |C,0  ]  =  U  +  y  1/2(1  +  C.  )U. 

i=2  1  1 

and 


K 

E[R-  |C,U  ]  =  l  1/2(1  -  C  )U  (5.30) 

i=2 


Since  the  desired  signal  has  an  average  count  'J^  the  average  counts  for  U  and 
■J  are 


K 

U  =  1/2  7  (1  +  C.)U. 
+  .-HO  i  - 


and 


K 

U  =  1/2  y  (1  -  C. )U. 

i=2  1  1 


(5.3D 
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To  uncondition  the  counting  probability  Pz(u  we  must  determine  the 

statistics  for  U+  and  U  .  This  is  an  intractable  problem  without  first  applying 

some  limiting  conditions.  For  example,  consider  the  expression  for  U+  which  is 

the  sum  of  K-1  independent  random  variables.  For  the  simple  case  of  a  single 

interferer  in  which  U.  is  gamma  distributed,  the  characteristic  function  for  the 

2. 

probability  density  of  U+,  conditioned  on  is 


M  (Ju>)  =  [1  -  jUl+Ci)<U.>/2mi  rmi 


(5.32) 


where  <IL>  and  m^  are  the  parameters  of  the  gamma  distribution.  For  N 
interferers,  the  characteristic  function  is 


N 

m  (ju)  s . n  [i 

+  1=1 


j'jj(1+  Cj.XUi>/2m  ]~mi 


(5.33) 


Mote  that  in  the  expression  for  M  (j'jj),  the  C^'s  are  all  independent  random 
variables  and  the  m  ' s  and  <  lh>'s  are  different  for  each  interferer.  As  a 
consequence,  it  is  difficult  to  invert  M  (ju))  to  obtain  the  probability  density 
function  of  U  .  The  sane  is  true  for  V  .  However,  we  can  simplify  the  problem 
considerably  if  we  assume  that  the  ' s  are  very  small,  as  is  the  case  for  long 
code  sequences,  and  the  fading  on  each  user  is  identically  distributed.  This 
latter  assumption  is  the  case  when  the  system  users  are  geographically  located 
in  a  region  characterized  by  similar  atmospheric  conditions.  Under  these 
conditions,  we  can  let  the  C^'s  =  0  and 
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<U^  =  <u2>  =  .  .  .  =  <uK>  =  <U> 

and 


In  this  case,  U+  is  equal  to  U  ,  so  that 


M+(joj)  =  M_(jw)  =  [1  -  joj<U>/2m]mN 


where  N  -  K-1 . 


The  probability  density  which  corresponds  to  the  characteristic 
given  in  C 5 . 34 )  is  the  gamma  density  with  parameters  N<U>/2  and  Nm 


p(U)  =  (2m/<U>)Nm  {UNm_1  /T (Nm) }  exp{-2mU/<U> } 


for  U  2  0, 

Under  these  simplifying  conditions,  (5.26)  becomes 


(5.34) 


function 


(5.35) 


(5.36) 
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Finally,  the  unconditioned  counting  probability  P^  is  determined  by 
averaging  (5-36)  over  the  gamma  probability  density  given  by  (5.35)*  The 
integration  gives  [Gradshteyn,  1965] 


P  (k)  =  F(  |k  ^Nm)/( j kj ) !r(Nm)  C  1  +  <U>/m  ]-tta 

•  (B/2)  k  2F1{(|k|+Nm)/2,  ( !  k  |+Nm+1  )/2;  |k|+1;32}  (5.37) 


where 

3  =  [  1  +  m/<U>  ]_1  (5.38) 


m  is  the  Nakagami  fading  parameter  and  F.|  is  a  hypergeometric  function  defined 
by 


2f1(  a,b;c;t 


=  T  {(a)  .(b)  /(c)  }  tj/j! 

j=0  J  J  J 


(5.39) 


where 


is  the  Pochhammer  symbol  for  the  product 


(a) j  =  a(a+1 )(a+2) . . . (a+j-1 ) ;  (a)Q  = 


1 


By  substituting  the  counting  probability  ?s  and  ?2  into  (5-21),  we  have 
shown  in  Section  D.2  of  Appendix  D  that  the  average  probability  of  error  equals 


P(E)  -  [1  ♦  <U>/m]-m(N+1)  I  C(Nm)r  /r! 3  (S/2)' 

r=0 


•  l  [(r+Nm)  /  p!]  ( @/2) p  I  /  s!'1  ^s.r-p 


(5 .  40) 


where 


N  is  the 


3  r  [  1  +  m/<U>  ]' 


number  of  interferers  and  m  is  the  fading  parameter.  <  3>r_p 


is  1/2  for 


s  =  r-p  and  is  1  otherwise. 


5.4  Downlink  Ghanne_l  Analysis. 


On  the  downlink  channel,  the  multiple  access  users  are  combined 
asychronously  and  coded  onto  a  single  output  signal;  therefore,  the  fading  is 
correlated  between  each  user.  Consequently,  the  average  probability  of  error 
for  users  on  the  downlink  can  be  determined  by  simply  deconditioning  the  error 
expression  obtained  for  the  no  fading  case  in  Chapter  4. 
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For  the  conditions  where  1)  the  cross-correlation  £  =  0,  2)  the  background 
radiation  and  detector  dark  currents  are  negligible  and  3)  the  fading  on  the 
intensity  is  slow  with  respect  to  the  data  bit,  the  error  expression  for  the 
downlink  (4.7),  generalized  to  a  conditional  probability  of  error,  is 


P(EIU)  =  1/2  [  1  -  Q(b’,a»)  +  Q(a',b')  ] 


(5.41) 


where  Q(a',b’)  is  the  Marcum's  Q-function  [Marcum,  I960],  and 


a’  =  [ (K ’  -  1)'J]1/2 


b'  =  t(K'  +  1 )U] 1 72 


(5.42) 


The  factor  K',  as  defined  in  (4.5),  is 


K'  =  K 


(  „ 


■  K/2 


(1/2) 


K 


K  even 


and 
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K-1 


K*  =  K 


(1/2)K-1 


|C  K— 1  ]/2J 


K  odd 


(5. 43) 


Notice  that  the  conditional  error  expression  for  the  downlink  is  identical 
to  that  for  correlated  fading  on  the  uplink,  with  K  replaced  by  K' .  Hence,  the 
average  probability  of  error  is  the  same  form  as  that  for  the  uplink  with  K 
replaced  by  K' 


P(E )  =  [1  +  K ' <U>/m]-m 


♦  l  (m)  [(K'-1 )/2K']k  ak  (k!)*1 

k=0 


•  l  (k+m),  [(K’+1)/2K']j  aj  (j!)-1  y,  ^  (5.44) 

j-0  J 


where 

a  =  [1  +  m/K,<U>]-1 

and  y.  =  1/2  for  j  =  k  and  is  1,  otherwise. 
J  t  K 
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5.5  Numerical  Results 

Error  probabilities  showing  the  effects  of  fading  on  optical  CDMA  systems 
are  presented  in  Figures  11-19  for  various  values  of  the  Nakagami  fading 
parameter  and  numbers  of  interfering  users.  The  results  are  computed  for 
negligble  background  noise,  =  0,  and  negligible  cross-correlation 
interference,  C.  =  0.  The  data  were  generated  by  numerically  evaluating  the 
error  expressions  for  correlated  (5.15)  an  uncorrelated  (5. *10)  fading  on  the 
uplink  and  for  correlated  fading  (5.1*1*)  on  the  downlink.  Unccrrelated  fading  is 
not  a  consideration  for  the  downlink  because  the  satellite  transmits  a  single 
optical  carrier  on  which  all  of  the  active,  uplink  users  have  been  modulated. 
Hence,  fading  of  the  optical  carrier  is  equivalent  to  each  user  undergoing  the 
same  instantaneous  fade;  therefore,  the  fading  is  correlated  between  each  user. 
On  the  other  hand,  correlated  fading  on  the  uplink  is  unrealistic  for  most 
systems;  since  this  implies  that  the  users  are  so  close  geographically  that  they 
all  experience  the  same  instantaneous  fades.  However,  this  case  was  easy  to 
solve  and  the  result,  with  a  simple  change  of  variable,  is  directly  applicable 
to  the  downlink.  In  all  of  the  numerical  evaluations,  series  expressions  were 
truncated  and  considered  a  satisfactory  approximation  after  summing  over  a 

_7 

sufficient  number  of  terms  to  insure  accuracy  to  within  10  .  For  greater 
accuracies,  the  uplink  expression  for  uncorrelated  fading  which  contains  the 
triple  summations,  converges  very  poorly;  especially  for  small  values  of  the 
fading  parameter  m  and  large  signal  intensities. 
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Figures  11 -13  examine  the  case  of  a  single  interferer  under  strong  fading 
conditions.  These  conditions  range  from  m  =  1 ,  the  Bose-Einstein  case  to  m  = 
10.  The  Bose-Einstein  case  characterizes  narrow  band  Gaussian  optical  fields. 
The  magnitude  of  the  optical  field  is  Rayleigh  distributed;  a  well-known  fading 
model  in  RF  systems.  The  error  rates  for  the  no  fading  case,  labeled  NF,  are 
included  in  each  figure  for  comparison.  As  m  increases,  the  error  rates  for  all 
fading  conditions  improve  as  expected.  The  uncorrelated  fading  on  the  uplink 


(Figure  12)  creates 

a  greater  degree  of 

degradation 

to 

the 

system  than 

correlated  fading  on 

either  the  uplink  or 

downlink . 

This 

is 

not  totally 

unexpected  because 

an  uncorrelated  fading 

condition 

here 

imp 

lies  tha"  the 

desired  signal  and  interferer  are  fading  independent  of  each  other.  The 
instantaneous  strength  of  the  fade  on  the  interferer  is  not  likely  to  be  the 
same  as  that  on  the  desired  signal.  This  randomness  is  an  additional  form  of 
noise  and  is  not  completely  averaged  out  in  tne  detection  process,  even  though 
the  cross-correlation  between  the  signals  is  negligible.  Whereas  in  the  case  of 
correlated  fading,  the  interferer  and  desired  signal  fade  instantaneously  and 
the  net  noise  effect  is  likely  to  be  less  after  the  averaging  in  the  detection 
process.  Over  the  performance  range  plotted,  the  downlink  exhibits  a  better 
error  rate  than  the  uplink  under  either  correlated  or  uncorrelated  fading 
conditions.  This  is  noted  by  the  wide  disparity  in  the  average  numbers  of 
photons  required  to  support  a  given  error  rate  for  a  specified  value  of  m. 
However,  the  system  as  a  whole,  suffers  serious  degradation  from  strong  fading. 
Much  larger  signal  intensities  are  required  to  obtain  reasonable  error  rates  for 


small  values  of  m. 


(SINGLE  INTERFERER) 


figure  IB.  Probability  of  error  vs.  average  signal  pliotons  for 
uncorrelated  fading  on  the  uplink  channel  for  a  single  interferer.  m  is 
the  Nakagami  fading  parameter.  Background  noise  and  cross-correlation 


'1 


(SINGLE  INTERFERER) 
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For  m  >>  1 ,  the  gamma  distribution  can  be  used  to  approximate  the  lognormal 
distribution.  Using  the  transformations 


ro  s  10  losiO(I) 

and 


(5.45) 


<iq>  =  io  iog1Q(<r>) 

Yoshikwa  (1979)  has  shown  that  the  probability  density  for  Iq  can  be 
approximated  by  the  form 


p(IQ)  =  l/vTT  - o  exp  C-(IQ 


<1  >)2  /2  32  ] 

0  o 


(5.46) 


where 


-  =  4.343/  m 

~o 

is  the  normalized  variance  in  terms  of  the  fading  parameter  m.  Hence,  with  I 
in  the  form  of  a  normal  distribution,  I,  the  intensity,  takes  the  form  of  a 


lognormal  distribution. 
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Figures  14-16  are  plots  showing  tne  effects  of‘  fading  on  the  multiple  user 
capabilities  of  the  uplink  and  downlink  channels,  for  cQ  -  .3  (m=210).  This 
fading  condition  is  characteristic  of  the  moderate  turbulence  measured  on  1000 
meter  terrestrial  links  by  Tatarski,  (1961)  and  Axford,  (  1979)  •  Figure  14  gives 
the  performance  of  the  uplink  channel  for  up  to  15  interferers,  under  correlated 
fading  conditions.  The  uncorrelated  f.  _  condition  is  given  in  Figure  15  for 
3  interferers.  Overall,  the  performance  of  the  uplink  seriously  degrades  when 
the  number  of  interferers  is  increased;  e  nccially  for  uncorrelated  fading.  For 
the  uplink  to  provide  reasonable  errcr  rates  under  fading  conditions  in  a 
multiuser  environment,  very  large  signal  intensities  are  required.  Comparing 
these  results  with  the  corresponding  error  rates  for  the  downlink  (Figure  16)  we 
note  that  the  downlink  error  rate  is  significantly  better  than  the  uplink  under 
these  moderate  fading  conditions,  as  a  much  larger  number  of  active  users  can  be 
supported  at  a  given  error  rate  and  photon  count.  This  is  advantageous,  because 
the  users  on  the  uplink  (ground  stations)  can  potentially  boost  their  signal 
levels  until  a  lower  error  rate  is  attained.  Whereas,  the  downlink  (satellite) 
is  power-limited.  The  additional  margin  of  performance  in  the  downlink  may  help 
offset  this  limitation. 

Figures  17-19  give  the  error  rate  performance  for  strong,  lognormal  fading, 
,7  =  .3.  Sy  comparison,  the  downlink  continues  to  perform  better  than  the 
uplink.  Overall,  increased  fading  causes  increased  degradation  to  the  system. 
Under  strong  fading  conditions,  such  as  the  lognormal  condition  in  Figures 
17-19,  it  is  unlikely  that  satisfactory  errcr  rates  from  end-to-end  can  be 
achieved  without  boosting  the  signal  intensities  on  the  uplink  to  unreasonable 
output  levels,  especially  to  support  a  large  number  of  users. 


Figure  VI.  Probability  of'  error  vs.  average  signal  photons  for 
correlated  fading  on  tin.-  uplink  channel  for  multiple  users  in  moderate 
lognormal  fading.  N  is  the  number  of  i  riter  ferers  and  o  is  the 
normalized  variance  of  the  1  >g- i  ritens  i  ty .  Background  noise  and 


CHAPTER  6 


CONCLUSIONS 


6 . 1  Conclusion 3 /Summary  of  Results 

In  this  study  we  have  investigated  the  effects  of  cross-correlation 
interference  and  intensity  fading  on  the  error  rate  performance  of  an  optical 
CDMA  communication  system.  The  analysis  is  restricted  to  direct  detection, 
optical  polarization  modulation  systems  employing  direct  sequence  modulation. 
We  have  developed  generalized  expressions  for  the  conditional  probability  of 
error  for  the  uplink  and  downlink  channels  of  a  multiaccess  satellite  repeater 
system.  The  expressions  are  given  in  terms  of  the  number  of  active  users, 
signal  strength  and  the  normalized  cross-correlation  between  user  codes. 
Numerical  results  were  obtained  for  systems  of  two  users,  for  both  length  31  and 
127  Gold  codes.  We  have  shown,  from  the  analysis  of  a  random  sampling  of  the 
length  3'  codes,  that  performance  is  significantly  affected  by  the  magnitude  of 
the  normalized  second  moment  of  the  cross-correlation  between  codes.  The 
smallest  error  rate  was  obtained  from  the  code  pair  having  the  smallest  second 
moment.  Thus,  in  a  multiaccess  system  in  which  the  code  length  is  fixed, 
interference  can  be  minimized  by  using  only  those  code  pair  combinations  which 
give  small  second  moments. 

We  have  also  shown,  in  comparing  the  performance  of  length  31  codes  to 
length  127  codes,  that  the  error  probability  decreases  with  increasing  code 
length.  This  is  not  unexpected;  for  as  the  code  length  increases  the  code 


properties  more  closely  approximate  those  of  the  purely  random  codes  for  which 
the  cross-correlation  between  codes  is  zero.  This  suggests  that  the  effects  of 
cross-correlation  interference  can  be  made  increasingly  small  by  using  longer 
length  codes. 

The  multiuser  capabilities  of  the  uplink  and  downlink  channels  were  shown 

to  differ  significantly.  Specifically,  it  was  shown  that  the  SNR  in  the  uplink 

is  asymptotically  proportional  to  1/K,  whereas  in  the  downlink  it  is 

1  /2 

asymptotically  proportional  to  1/K  .  This  means  that  for  K  users  at  a  given 

2 

error  rate  in  the  uplink,  K  users  can  be  supported  at  the  same  error  rate  in 

the  downlink.  Therefore,  the  overall  system  capacity  is  limited  by  the  capacity 

of  the  uplink.  However,  this  problem  can  be  solved,  to  some  extent,  by  dividing 

the  uplink  coverage  area  into  geographically  separated  cells.  These  cells  could 

be  spatially  multiplexed  in  the  satellite  receiver  by  imaging  each  cell  onto  a 

photodiode  element  of  an  array  of  photodiodes.  Then,  by  multiplexing  K  cells  of 

2 

K  users,  we  could  extend  the  capacity  of  the  uplink  to  K  users.  This 
effectively  matches  the  capacity  of  the  downlink  for  a  comparable  error  rate. 
However,  in  this  case  "dead  spots"  between  adjacent  elements  of  the  photodiode 
array  would  introduce  some  degradation. 

In  the  analysis  of  fading  effects  on  optical  CDMA  systems,  we  have  used  the 
Nakagami  m-distribution  as  a  generalized  statistical  model  for  fading  on  the 
received  optical  intensity.  The  average  probability  of  error  was  derived  for 
both  correlated  and  uncorrelated  fading  in  terms  of  the  number  of  interfering 
users ,  the  Nakagami  fading  parameter  and  the  average  number  of  photons  counted 
in  a  bit  interval.  The  results  are  valid  for  a  noise-free  system  in  which  the 
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cross-correlation  interference  is  negligible.  Numerical  results  were  obtained 
for  both  strong  and  moderate  lognormal  fading,  using  the  Nakagami  m-distribution 
to  approximate  the  lognormal  distribution.  The  results  show  a  significant 
difference  between  the  multiuser  capabilities  of  the  uplink  and  downlink 


channels  under  fading  conditions 

.  It  was  shown 

that 

the 

downlink 

performs 

consistently  better  than 

the 

uplink  under 

all 

of 

the 

fading  conditions 

evaluated.  This  margin  of 

performance  derives 

from 

the 

SNR 

advantage 

of  the 

downlink,  i.e.,  the  SNR  of  the  downlink  is  asymptotically  proportional  to 
1  /2 

1/K  while  in  the  uplink  is  asymptotically  proportional  to  1/K.  As  a 
consequence,  the  overall  system  capacity,  under  fading  conditions,  is  limited  by 
the  capacity  of  the  uplink. 

Strong  fading,  which  is  characterized  by  small  values  of  the  fading 
parameter  (m  <  5),  seriously  degrades  the  performance  of  this  system.  However, 
moderate  lognormal  fading  with  intensity  variances  comparable  to  the  measured 
data  in  the  references  [Tatarski,  1961]  and  [Axford,  1979],  introduces  moderate 
degradation  and  reasonable  error  rates  for  multiple  users  can  be  obtained  by 
boosting  the  signal  power. 

6 . 2  Recommendations  for  Future  Study 

The  general  study  of  cross-correlation  interference  in  optical  CDMA  systems 
offers  many  logical  extensions.  Some  which  closely  parallel  this  study  are, 
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1)  Development  of  methods  to  improve  the  multiuser  capability  of  the  uplink 
channel . 

The  foregoing  analysis  showed  that  the  weak  link  in  the  end-to-end 
performance  of  this  system  is  the  uplink  channel.  In  spite  of  setting  the 
cross-correlation  to  zero,  each  additional  user  to  the  uplink  contributed  a  shot 
noise  count  equal  to  U/2.  This  noise  count  is  directly  proportional  to  the  chip 
period  of  the  multiple  access  code  sequence.  Therefore,  an  approach  to  reducing 
this  noise  contribution  is  to  transmit  shorter  pulses;  that  is,  instead  of 
transmitting  a  constant  carrier  for  the  full  chip  period,  transmit  over  only  a 
fraction  of  the  period.  The  receiver,  gated  to  correlate  over  only  a  fraction 
of  the  chip  period,  would  less  likely  be  gated  on  at  the  same  time  that  the 
interfering  signal  is  present  at  the  receiver  input.  We  expect  that  if  the 
carrier  is  gated  to  1/100-th  of  a  chip  period,  for  example,  the  net  effect  would 
be  a  1/100-th  reduction  in  the  shot  noise  contribution. 

2)  Development  of  a  simple,  easy  to  evaluate,  approximation  to  the 
probability  of  error  in  terms  of  the  moments  for  the  cross-correlation  values. 

From  a  system's  point  of  view,  it  would  be  invaluable  to  have  a  simple 
expression  relating  the  gross  behavior  of  the  error  probability  as  a  function  of 
the  number  of  users,  moments  for  the  cross-correlation  values  and  the  signal 
intensity.  System  designers  could  then  obtain  an  order  of  magnitude  feel  for 
the  code  lengths  and  signal  powers  required  to  support  a  given  number  of  users 
at  a  specified  error  rate. 
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3)  Analysis  of  the  error  rate  performance  of  optical  heterodyne  detection 
systems . 

Optical  heterodyning  is  accomplished  by  adding  the  received  intensity  to  a 
strong  locally  generated  reference  intensity  and  photodetecting  the  sum.  The 
prime  objective,  as  with  conventional  heterodyning,  is  to  generate  a 
frequency-translated  version  of  the  received  field  at  a  frequency  which 
optimizes  detection.  If  the  local  reference  is  strong  enough,  photoelectrons 
are  emitted  from  the  detectors  at  rates  which  are  so  high  that  the  detector 
output  signal  is,  in  effect,  a  replica  of  the  received  field  plus  an  additive 
Gaussian  noise  contribution  from  the  reference  intensity.  In  this  case,  the 
primary  source  of  detection  noi3e  is  that  induced  by  the  reference  intensity; 
the  shot  noise  generated  by  the  signal  and  detector  dark  currents  can  be 
neglected.  SNR  comparisons  of  heterodyne  and  direct  detection  systems  show  that 
for  certain  frequencv  ranges  heterodyning  is  preferred  to  direct  detection 
[Fried,  1967].  Extending  this  study  to  include  an  analysis  of  optical 
heterodyne  systems  would  provide  a  performance  comparison  between  direct 
detection  and  heterodyne  CDMA  systems  in  the  presence  of  cross-correlation 
interference  and  form  a  basis  for  choosing  the  appropiate  receiver  for  a 
specific  application. 

4)  Extension  of  the  analysis  technique  to  the  study  of  cross-correlation 
interference  in  optical  ranging  systems  which  use  direct  sequence  modulation. 
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Range  measurements  are  made  by  determining  the  time  required  for  a  signal 
to  propagate  to  an  object  and  return.  Optical  ranging  using  DS  modulation  can 
be  accomplished  by  modulating  the  polarization  of  the  transmitted  signal  with  a 
pseudorandom  code  sequence  (ranging  code) .  The  return  signal  is  a  time-delayed 
version  of  the  transmitted  signal  plus  noise.  The  receiver,  using  correlation 
techniques,  extracts  the  relative  time  delay  and,  indirectly,  the  range 
information.  The  analysis  procedure  applied  to  cross-correlation  interference 
is  also  applicable  to  an  optical  DS  ranging  system.  However,  in  this  case,  the 
autocorrelation  properties  of  a  particular  ranging  code,  as  opposed  to  the 
cross-correlation  between  different  codes,  is  an  important  system  consideration. 
By  appropiately  modeling  the  photoelectron  emission  rate  of  the  return  signal, 
we  can  analyze  the  effects  of  atmospheric  fading,  doppler  shift,  signal 
depolarization,  or  background  radiation  on  range  resolution  and  accuracy. 
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APPENDIX  A 


COUNTING  STATISTICS  FOR  THE  RECEIVER  OUTPUT 

In  this  appendix  we  derive  an  alternative  expression  for  the  probability 
distribution  of  the  photon  count  at  the  receiver  output  (3*4)  in  Chapter  3*  The 
receiver  output  is  analyzed  as  the  diferenoe  between  two  independent  Poisson 
random  variables. 

Let  R  =  R+-  R~  be  the  difference  between  two  independent  Poisson  random 
variables  (  R+  and  R“) .  Then  the  probability  distribution  of  R  is  given  by 
Pratt  (1969) 

PR(R=n)  =  exp[-(E[R+3+E[R'*] )  {E[R"]/E[R+] }"n/2 

I|n| {2  E[R+3  E[R']}  (A1) 

where  E[  •  ]  denotes  expectation.  To  express  the  distribution  of  R  in  terms  of 
its  mean  and  variance,  we  note  that 


E[R]  =  E[R+3  -  E[R”3 


( A2) 


and 


(A3) 


Var(R)  s  Var(R+)  +  Var(R") 
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where  we  have  used  the  fact  that  the  Poisson  random  variables  are  independent. 
Since  the  mean  and  variance  of  Poisson  random  variables  are  equal,  (A3)  can  be 
written: 


Var(R)  =  E[R+]  +  E[R“] 


(A4) 


Solving  (A2)  and  (A4)  for  E[R+]  and  E[R~],  we  get 


E[R+]  =  1/2  (  Var(R)  +  E[R]  } 


(A5) 


ElR-]  =  1/2  {  Var(R)  -  E[R]  }  (A6) 

The  probability  distribution  of  R  can  now  be  written: 


PR(R=n)=  exp  { -Var(R)} { (Var(R)  -  E[R] )/( Var(R)  +  E[R])}'n/2 

1  |n,  (tVar(R)2  -  E[R]2]1/2}  (AT) 
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APPENDIX  B 

DERIVATION  OF  THE  CONDITIONAL  PROBABILITY  OF  ERROR 

In  this  appendix  we  derive  the  conditional  probability  of  error  for  an 
arbitrary  number  of  independent  users  with  different  signal  intensities.  The 
error  expression  is  conditioned  on  the  cross-correlation  C^,  which  is  defined  as 
the  cross-correlation  between  the  code  of  the  desired  user  and  the 
code-modulated  data  signal  of  the  i-th  interferer.  This  derivation  is 
facilitated  by  the  fact  that  the  counting  statistics  of  the  receiver  output  can 
be  expressed  in  terms  of  its  mean  and  variance.  This  is  shown  in  Appendix  A. 
Therefore,  as  a  first  step,  we  derive  the  mean  and  variance  of  the  receiver 
output,  defined  by  (3*0.  We  designate  the  desired  signal  as  b^t)  with  code 
sequence  a^t),  and  to  simplify  notation,  use  the  designation  E[R]  to  represent 
the  conditional  mean  ECRIC^. 

B. 1  Mean  and  Variance/Uplink  Channel 

For  K  independent  users,  the  photodetector  outputs  can  be  written 

K 

Y(t)  =  l  Y.  ( t)  (B1 ) 

i-1  1 

and 

K 

X(t)  =  l  X.(t)  (B2) 

i-1  1 
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where  Y^t)  and  X^(t)  are  the  photodetector  emissions  due  to  the  i-th  user.  The 
time-varying  rate  parameters  of  the  emission  processes  are  proportional  to  the 
instantaneous  polarity  of  the  coded  signals.  These  rate  parameters  are 


fXiai(t-xi)bi(t-Ti)  +Xfe/K  a^  >  0 


Xfe/K  aibi  <0 


X. /K  a.b,  >  0 

b  i  J 

-Xia1(t-Ti)bi(t-xi)  +  Xb/K  aibI  <  0  (B3) 

where  X  ^  is  the  photon  counting  rate  of  the  i-th  user  and  X  b  is  the  rate 
parameter  for  the  background  noise,  shown  equally  distributed  among  the  K  users. 
The  conditional  mean  of  (3.1)  is 


E[R] 


a1(t){E[Yi]  -  E{X1]}dt 


(B4) 


Since  the  mean  value  of  a  Poisson  process  is  equal  to  the  rate  parameter,  we 


have 
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EtR]  =1 
i=l 


K  ,Td 


a.(t)[X  -  X  ]dt 

n  1  ” .  x . 

u  11 


Substituting  (B3)  into  (B5),  we  get 


E[R] 


K  r  d 

=  I  Xi  a, 

i=l  iJ  0  1 


(t)ai(t-xi)  bi(t-xi)dt 


2 

This  is  simplified  by  noting  that  a*(t)  is  identically  one  and  b1  is 
over  the  integration  period.  The  conditional  mean  reduces  to 


K 


E[R]  =  b  U  +  I  U.C. 

i  1  £=2  1  1 


where 


,Td 


Ci  =  1/Td 


a1(t)ai(t-r1)  bt(t-ri)dt 


is  the  normalized  cross  correlaton  and  U.  =  X.T. . 

1  la 


The  variance  is  obtained  in  a  similar  manner 


K  K  rTdrTd 


Var(R)  =  l  l 

j=l  i=l 


a1(t1)a1(t2)  {ECYi(t1)YJ(t2)3 


+  E[Xi(t1)X  (t2)]  -  E[Xi(t1)Yj(t2)] 


-  E[Xj(t2)T1(t1>])  dt1dt2 


(B5) 


(B6) 

constant 


(B7) 


(B8) 


(B9) 
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where  for  convenience  we  have  defined 


Yt(t)  =  Yi(t)  -  E[Yi(t)] 


and 


X1(t)  =  Xj,(t)  -  E[X  (t)] 


Under  the  assumption  that  the  user  signals  are  mutually  independent 
independent  of  the  background  noise,  (B9)  simplifies  to 


K  fTd  fTd 


rv  u  f  a 

Var(R)  =  l  a  (t'  )a  (t  )[  C  (t.,  t_) 

i=]_J  Q  J  Q  '  '  1  2  y  1  2 


cx  <t,, 


Cy(t^,  t^)  and  C^Ct^,  are  the  covariance  functions  associated 
Poisson  processes  Y^t)  and  Xt(t).  They  are  given  by  [Papoulis,  1965] 


c(tr  t2)  =  X(t2)fi(t1  -  t  ) 


Substituting  (B12)  into  (B11),  we  get 


Var(R) 


K 

■Ji 


Ld 

[\ 

0  yi 


V]dt 

l 


(BIO) 

and  also 

(B11) 

with  the 

(B12) 

(B13) 


From  (B3)  we  obtain 
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X  +Xx  =  2Xb  /K  +  Xj_  |ai(t-xi)bi(t-Ti)|  (B14) 

i  i 

The  factor  |  |  is  identically  one  for  all  values  of  a.^; 

therefore,  the  variance  reduces  to 


K 


Var(R)  =  2U.  +  y  U. 

b  i=l  1 


where  U.  =  X.T.  and  U.  =  X.T. . 
i  id  b  b  d 


B.2  Mean  and  Variance/Downlink  Channel 


(B15) 


For  the  downlink  channel,  the  rates  at  which  photoelectrons  are  emitted 
from  the  photodetectors  in  the  receiver  are  proportional  to  the  instantaneous 
magnitude  of  the  multiplexed  signals.  These  rates  are 


K 

xl  »1(t-T i)b i(t-xi)  +  Xfe 

i=l 


l  aibi  >  0 

y  a. b.  <  0 

L  li 


aibi 


>  0 


X 


x 


K 

-X  l  a  (t-T.)b.(t-T  )  +  Xfe  l  a  b  <  0  (B16) 

>  1=1 
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where  A  is  a  constant  photon  counting  rate  which  is  multiplied  by  the 
instantaneous  magnitude  of  the  multiplexed  signals  to  form  the  magnitude  of  the 
transmitted  signal. 

The  conditional  mean  of  R  is  determined  in  a  manner  similar  to  that  of  the 
previous  case.  The  result  is 


K 

BCR]  =  b.U  +  U  l  C.  (B17) 

1  1=2  1 


where  U  =  A  T . . 

a 


If  we  apply  the  same  independence  assumptions  that  were  used  in  the 
previous  case,  the  variance  reduces  to 


Var(R) 


[A  (t)  +  A  (t)]dt 

0  y  x 


(B18) 


Substituting  ( B 1 6 )  into  (B18),  we  get 


or 


Var(R) 


T 

d  K 

[2A  +  A  |  £  a.  (t-t. )  b.(t-t,  )j]dt 

o  _  1  i  1  1  i 

0  1  -1 


( B 1 9 ) 


Var(R)  =  2U.  +  K'U 
b 


( B20 ) 
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where 


fTd  K 

K’  a  1/T.  |  l  a.(t-T.)b.(t-T,)  |dt  (B21) 

dJ0  i=l  1  11  1 


B.3  Derivation  of  the  Error  Probability 

If  b1  a  +1  are  equally  likely,  the  conditional  probability  of  error, 
defined  at  the  receiver  output  R  is 


P(E  C.)  =  1/2Pr(R>0j  Cit-1)  +  1/2Pr(R<0|CiI  +  1) 

1/4Pr(R=0  |C± ,-1 )  +  1/4Pr(R=0  |C  ,+1)  (B22) 

where  is  the  cross-correlation  given  in  (B8).  Denote  the  conditional 
counting  probability  of  the  receiver  output  by  P  (k  |c. ,b.),  then  the  conditional 

n  11 

probability  of  error  can  be  written 


00  00 

P(E |Ci)  =  1/2[1  ♦  l  PR(k|Clt  -1)  -  l  PROe|C1,  +1) 

0  1c*5  0 

+  1/2PR(o!ci(  +  1)  -  1/2  PR (0 1 Ci ,—!)]  (B23) 


From  Appendix  A,  P_  can  be  expressed 

n 
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PRCk  |Ci)  =  exp  {-Var(R)}  {(Var(R)-ErR])/(Var(R)+E[R])}_k/2 


I  |  k  j  {/var(R)2  -  E[R]2  } 


(B24) 


Consider  first  the  derivation  of  the  uplink  probability  of  error.  By 
substituting  the  mean  (B7)  and  the  variance  (B15)  of  R  into  (B24),  we  obtain  the 
receiver  counting  statistics  for  the  uplink 


P(k|ct,  +  1)  =  exp{-(a2  +  b2)/2}  (b/a)k  X  |k(Cab) 


(B25) 


and 


P(k  C± ,-1 )  =  exp{-(c2  +  d2)/2}  (d/c)k  I^Ccd) 

where 


(B26) 


a  =  [  2Ub 
b  =  [  2Ub 
c  =  [  2Ub 
d  =  [  2Ub 


K  1/? 
+  Y  U.(1-C.)]1/2 

i=2  1  1 


2U. 


2U. 


K 

+  Y  U. ( 1+C. ) ] 
i=2  1  1 


K 

+  I  u.d-c,)] 

i=2  1  1 


1/2 


1/2 


*  l  uid+ci)]1/2 
i=2 


(B27) 
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When  (B25)  and  (B26)  are  substituted  into  (B23)  and  we  use  the  series  expansion 
of  the  Marcum's  Q-function  [Marcum,  I960]  given  by 


Q(a,b)  =  exp{-(a2  +  b2)/2)  I  (a/b)j  I.(ab) 

j=0  J 


to  simplify,  we  get 


(B28) 


P(EIC1)  =  1/2C1  -  Q(b,a)  +  Q(a,b)] 


K 


+1 


/4exp{-(2U  *  l  U,)HI  (ab)  -  I  (cd)} 

D  ^  i  o  O 


(B29) 


where  the  arguments  of  the  Q-function  are  given  in  (B27). 


The  probability  of  error  for  the  downlink  is  derived  in  a  similar  manner 
but  using  the  mean  (B17)  and  the  variance  (B20)  of  the  receiver  output  for  the 
downlink.  The  significant  difference  however,  is  that  we  must  account  for  the 
fact  that  K'  is  also  conditioned  on  b^(t)  and  C^.  Hence,  in  the  generalized 
expressions  for  the  conditional  probability  of  error  for  the  downlink  we  denote 
this  conditioning  by  defining 


and 


K'  *  given  b.  s  +1 
+  I 


K '  s  K'  given  b  =  -1 


This  dependence  is  shown  explicitly  in  the  evaluation  of  K'  for  some  specia 
cases  in  Appendix  C.  The  conditional  probability  of  error  for  the  downlink  is 


P(E|Cl)  =  1/2[1  -  Q(b'.a')  +  Q(d',c')] 

+1/4exp[-(2U.  +  K'U)]  I  (a'b*) 

D  +  0 

-1/4exp[-(2U.  +  K'U]  I  (c’d')  (B30) 

D  -  O 


where 
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APPENDIX  C 


SIMPLIFICATION  OF  THE  EXPRESSION  FOR  K' 


In  this  appendix,  we  derive  two  simplified  expressions  for  K':  1)  for  a 
single  interferer  and  2)  for  an  arbitrary  number  of  interferers  when  the 
cross-correlation  equals  zero. 

c.i  single  latficiffican 


For  a  single  interferer,  K'  is  given  by 


T. 

r  <3 


1/T, 


a1  (t) 


b1  (t) 


♦  a2(t-r2) 


b2(t-t2)  |dt 


(Cl) 


where,  since  time  delays  are  referenced  to  a^t),  we  have  sett  ^x  0.  If  a1b1 
and  a2b2  are  of  the  same  polarity 


a1  b1  +  a2  b2  =  2  (C2) 

A  A 

On  the  other  hand,  if  and  *^2  are  °PP°site  Polarity 


al  b1  ♦  a2  b2  *  0  (C3) 


We  can  conveniently  express  these  relations  with  the  single  form 
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UaiV2b2 


(CM) 


Substituting  the  above  into  (Cl),  we  get 


K' 


1  /T  . 


(1  ♦  a1(t)b1(t)  a2(t-T2)b2(t-x2)]dt 


(C5) 


Since  b^(t)  is  constant  over  the  bit  interval,  it  can  be  brought  outside  the 
integral.  Performing  the  integration  gives 


,Td 


a^  (t) 


b1(t)a2(t-x2)b2(t-^)dt] 


(C6) 


The  factor  in  braclcets  is  the  normalized  cross-correlation  C2  between  a^  and 
a2b2  for  delay  t2-  Therefore,  K'  takes  the  particularly  simple  form 


K'  *  1  ♦  b,C2 


(C7) 


c.2  Muiualc  tot&rjfcr.gf.a 

In  the  limit  as  C,  ♦  0,  the  general  expression  for  K'  takes  the  form  of  a 
time  average  of  the  absolute  value  of  the  sum  of  finitely  many  independent 
random  variables  which  take  on  the  values  +1  and  -1.  The  chance  fluctuations  in 
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the  value  of  the  summation  in  the  integrand  of  JO  can  be  modelled  analytically 
as  a  random  walk  process  with  the  integrand  taken  to  be  the  absolute  value  of 
the  path  length  from  the  origin  to  the  final  position  of  a  man  or  particle 
performing  a  random  walk  [Feller,  1966].  Assuming  a  statistically  stationary 
system,  JO  can  then  be  analyzed  as  the  expectation  of  the  absolute  value  of  the 
path  length  after  K  random  steps,  where  K  corresponds  to  the  number  of 
simultaneous  users  in  the  downlink  channel. 

Let  be  the  path  length  after  K  steps  then  the  probability  distribution 
of  |  Z  )  is  given  by 


Pr(  \Zk\  «n) 


(JUn)/2 


(1/2)' 


CC8) 


where  n  a  0,1,2,...,  K,  and  the  binomial  coefficient  is  taken  to  be  zero  unless 
(K+n)/2  is  an  integer  <  K.  The  Neumann  factor  e  ,  is  1  if  n=0  and  2 


otherwise. 
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For  K  even  we  have 


K'  *  E[(  zj 


K/2 

:  [  k 
n-1 


Un+K/2)  J 


(1/2)K 


(C9) 


which  by  a  simple  change  of  variable  can  be  written  as 


K 

K 

K 

K 

c 

J 

ja¬ 

il 

(1/2)K  -  2K  T 

(1/2)K 

n-K/2+l 

.n  ^ 

n-K/2+1 

I 

A 

(CIO) 

Using  the  properties  of  the  binomial  distribution  to  simplify  the  summations 
the  above  equation,  we  get 


K'  *  K 


K/2 


(1/2)' 


CC11) 


in 


.  > 
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The  derivation  for  K  odd  follows  in  a  similar  manner;  the  result  is 


r 


K'  *  K 


K-1 


(i/a)*-1 


(.  (K-1  )/2j 


(C12) 


APPENDIX  D 


THE  PROBABILITY  OF  ERROR  FOR  CORRELATED  AND  UNCORRELATED  FADING 


In  this  appendix,  we  derive  simplified  expressions  for  the  average 
probability  of  error  for  multiple  users  on  the  uplink  channel,  for  correlated 
and  uncorrelated  fading. 

D. 1  Correlated  Fading 

The  average  probability  of  error  for  correlated  fading  on  the  uplink 
channel  is  obtained  by  by  integrating  the  conditional  probability  of  error 
(3.18)  over  the  probability  density  for  the  fading  intensity 


P(E) 


P(EIU)  p(U)  dU 


(D 1 ) 


where  p(U)  is  the  probability  density  for  the  fading  intensity,  given  by 


p(U)  =  (m/<U>)m{Um“1  /T(m)}  exp{-m(U/<U>) } 


(D2) 


for  U  >  0  and 


P<E|U)  »  1/2  [  I  -  Q(b,a)  +  Q(a,b)  ] 


(D3) 


where  Q(a,b)  is  the  Marcum's  Q-function  [Marcum,  I960],  with  arguments 


Ill 


a=  [(K-  1)U]1/2 


and 


b=  [(K+  1)0]1/2 


(D4) 


Using  the  symmetry  relation  [Helstrom,  I960] 


Q(a,b)  +  Q(b,a)  =  1  +  exp{-  (  a2  +  b2)/2)  IQ(ab)  (D5) 


the  series  expansions  for 


Q(a,b)  =  exp{-  (  a2  ♦  b2)/2}  J  (a2/2)k  /k!  \  (b2/2)J  /J! 

k«0  j-0 

(D6) 


and  the  modified  Bessel  function 


I.(a)  =  7  (a/2)2k  /k!k!  CD7) 

0  k=0 


in  (D3).  we  get 


00 

P(E  U)  =  exp(-KU)  J  [U(K  -  1)/2]P/p! 
p=0 

•  l  [U(K  +  1)/2]J/J«  Y , 
j-0  J,P 


(D8) 
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where  Y .  =  1/2  for  jsp  and  1  otherwise. 

J  ♦  P 


If  (D8)  and  (D2)  are  substituted  into  (D1)  the  required  integration 


°o  p 

P(E)=(m/<U>)m/r (m)  l  l  [(K  -  1)/2]P/p  C(K  +  1)/2]j/j!  Y 
p=0  j=0 


*  [  uJ"‘P+m"1  exp{- (K  +  m/<U>)U}  dU 


Using  the  integral  identity  [Gradshteyn,  1965] 


.00 

xP_1  exp(-qx)  dx  s F  (p)  q-p 

J0 


we  get 


P(E)  =  [1  +  K<U>/m  ]’m 


•l  [(K  -  1)/2K]P6  p  [(m)  /p! ] 
p=0  p 


.  V  [  (K  +  1)/2K]j  5J  [(m+p)  ,/j!  ]  y . 

j  =0  J  J  >P 


where 


<5  =  [  1  +  m/X<U>  ] 


becomes 


(D9) 


(DIO) 


(Dll) 


(D12) 
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D.2  Uncorrelated  Fading 

The  equation  for  the  average  probability  of  error,  expressed  in  terms  of 
the  counting  probabilities  for  the  desired  signal  S  and  the  interference  Z,  is 
(5.20) 


(D13) 


where,  y.  .  is  1/2  for  j=k  and  1  otherwise.  P  is  the  desired  signal  counting 

j  s 

probability  (5.24) 


Ps(j)  =  F (m+j)/F(m)  j!  [1  +  <U>/m  ]"m  SJ 

for  j  2  =  0.  (D 14 ) 

where  <'J>  is  the  mean  signal  count,  m  is  the  Nagakami  fading  parameter  and 


8  =  [  1  +  m/<U>  l"1  ( D 1 5 ) 

Since  Pg  is  zero  for  J  <  0,  the  summation  is  (D13)  is  reduced  to  a  summation 

over  the  positive  range  of  the  index  k.  P  is  the  counting  probability  for  the 

z 

interference,  which  for  the  positive  integers  is  (5.31) 
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P  (k)  =  r (k+Nm) / r(Nm)  k!  [1  +  <U>/m]"Nm  (B/2)k 
z 

•  2Ft{(k+Nm)/2,  (k+Nm+1 )/2;  k+1;  S2}  (D16) 

where  N  is  the  number  of  interferers 


6=  [  1  +  m/<U>  ) 

and2F1  is  a  hypergeometric  function  defined  by  Abramowitz,  (1965) 


F  (a,b;c;t)  =  £  {(a) .(b)  /(c)  }  tp/p! 

*  1  p=0  p  p  p 


(D17) 


The  notation  (a)p  is  the  Pochhammer  symbol  for  the  product 


(a)  =  a(a+1 ) (a+2) . . . (a+p-1 ) ;  (a)  =  1 

P  0 

The  Nakagami  fading  parameter  m  and  the  mean  signal  count  <U>  are  the  same  as 
those  for  the  desired  signal. 


By  substituting  the  counting  probabilities,  ?s  and  Pz>  into  the  error 
expression,  we  get 


P(E)  =  [  1  +  <U>/m  ] 


-m(N+1 ) 


l  (Nm)  /p!(B/2)! 


p=0 


?  P  1 

•  ,F1{(p+Nm)/2,  (p+Nm+1 )/2;  p+1;  B  }  l  (m)./j!  6 J  y,  „ 

£  j=0  J  J,p 

(D1S) 
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where  we  have  used  the  identity  [Hansen,  1975] 


(m)^  =  r(m+k)/T(m) 


to  simplify  notation. 


Using  the  relation  [Abramowitz,  1965] 

iT1/2r(2p)  =  22p_1  r(p)r(p+i/2) 


we  can  write  the  hypergeometric  function  in  ( D 1 8 )  as  follows, 


2F1{(p+Nm)/2,  (p+Nm+1)/2;  p+1;  32}  = 


1 / IXp+Nm)  l  [F (p+Nm+2r)/(p+1)  r!  ]  (B/2)2r 

r=0  r 


Substituting  (D21)  into  (D18)  and  applying  the  identity  [Hansen,  1975] 


00  00  00  p 

I  I  f(r,p)  =  l  l  f(r,p-r) 

p=0  r=0  p=0  r=0 


to  reorder  the  summations,  we  obtain 


PCS)  =  [  1  +  <U>/m  ] 


•  3r+P+j  (  i/2) 


-■(N+1 ) y  f  ryP 

r£0  p£0  j£0  7J.r-p 

r+P  [f  (p+Nm+r)/r(Nm)]  [  (m)  / j ! p!  r !  ] 


(D19) 


(D20) 


(D21 ) 


(D22) 


(D23) 
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This  can  be  reduced  further  using  the  identity  ( D 1 9 ) ;  the  resulting  average 
probability  of  error  is 

00 

P(E)  r  [1  +  <U>/m]_(N+1)m  l  (Mm)  /rl  (g/2)r 

r=0  r 

•f  (r+Nm)p/p!  (g/2)pr^P(m)  j/j !  6jYj>r_p  (D2«) 
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